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Abstract 



Analogue models of (and for) gravity have a long and distinguished 
history dating back to the earliest years of general relativity. In this 
review article we will discuss the history, aims, results, and future 
prospects for the various analogue models. We start the discussion by 
presenting a particularly simple example of an analogue model, before 
exploring the rich history and complex tapestry of models discussed 
in the literature. The last decade in particular has seen a remarkable 
and sustained development of analogue gravity ideas, leading to some 
hundreds of published articles, a workshop, two books, and this review 
article. Future prospects for the analogue gravity programme also look 
promising, both on the experimental front (where technology is rapidly 
advancing) and on the theoretical front (where variants of analogue 
models can be used as a springboard for radical attacks on the problem 
of quantum gravity). 
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1 Introduction 



And I cherish more than anything else the Analogies, my most 
trustworthy masters. They know all the secrets of Nature, and 
they ought least to be neglected in Geometry 

— Johannes Kepler 

Analogies have played a very important role in physics and mathemat- 
ics - - they provide new ways of looking at problems that permit cross- 
fertilization of ideas among different branches of science. A carefully chosen 
analogy can be extremely useful in focussing attention on a specific problem, 
and in suggesting unexpected routes to a possible solution. In this review 
article we will focus on "analogue gravity", the development of analogies 
(typically but not always based on condensed matter physics) to probe as- 
pects of the physics of curved spacetime — and in particular to probe aspects 
of curved space quantum field theory. 

The most well-known of these analogies is the use of sound waves in a 
moving fluid as an analogue for light waves in a curved spacetime. Super- 
sonic fluid flow can then generate a "dumb hole", the acoustic analogue of 
a "black hole", and the analogy can be extended all the way to mathemat- 
ically demonstrating the presence of phononic Hawking radiation from the 
acoustic horizon. This particular provides (at least in principle) a concrete 
laboratory model for curved-space quantum field theory in a realm that is 
technologically accessible to experiment. 

There are many other "analogue models" that may be useful for this or 
other reasons — some of the analogue models are interesting for experimental 
reasons, others are useful for they way they provide new light on perplexing 
theoretical questions. The information flow is in principle bi-directional and 
sometimes insights developed within the context of general relativity can be 
used to understand aspects of the analogue model. 

Of course analogy is not identity, and we are in no way claiming that the 
analogue models we consider are completely equivalent to general relativity 
- merely that the analogue model (in order to be interesting) should capture 
and accurately reflect a sufficient number of important features of general 
relativity (or sometimes special relativity). The list of analogue models is 
extensive, and in this review we will seek to do justice both to the key 
models, and to the key features of those models. 



S 



The plan for this review is as follows. In the following chapters we shall: 

• Discuss the flowing fluid analogy in some detail. 

• Summarize the history and motivation for various analogue models. 

• Discuss the many physics issues various researchers have addressed. 

• Provide a catalogue (hopefully complete) of extant models. 

• Discuss the main physics results obtained to date. 

• Outline the many possible directions for future research. 

• Summarize the current state of affairs. 

By that stage the interested reader will have had a quite thorough in- 
troduction to the ideas, techniques, and hopes of the analogue gravity pro- 
gramme. 

1.1 Going further 

Apart from this present review article, and the references contained herein, 
there are several key items that stand out as starting points for any deeper 
investigation: 

• The book "Artificial Black Holes", edited by Mario Novello, Matt 



Visser, and Grigori Volovik [|286| . 

The websites for the 'Analogue models" workshop: 

— http: / /www. cbpf.br/~bscg/analog/ 

— http:/ /www.mcs.vuw.ac.nz/~visser/Analog/ 

— http:/ /www. physics. wustl.edu/~visser/Analog/ 



The book "The Universe in a Helium droplet" , by Grigori Volovik |419| . 



The Physics Reports article, "Superfluid analogies of cosmological phe- 
nomena" , by Grigori Volovik [|414| . 
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2 The simplest example of an analogue model 

Acoustics in a moving fluid is the simplest and cleanest example of an ana- 
logue model | |388| , [390| , |391|| . The basic physics is simple, the conceptual 



framework is simple, and specific computations are often simple (whenever, 
that is, they are not impossibly hard). 1 

2.1 Background 

The basic physics is this: A moving fluid will drag sound waves along with it, 
and if the speed of the fluid ever becomes supersonic, then in the supersonic 
sound waves will never be able to fight their way back upstream ||378| , [388| , 



i90[ |391| | . This implies the existence of a "dumb hole" , a region from which 
sound can not escape. 2 Of course this sounds very similar, at the level of a 
non-mathematical verbal analogy, to the notion of a "black hole" in general 
relativity. The real question is whether this verbal analogy can be turned 
into a precise mathematical and physical statement — it is only after we have 
a precise mathematical and physical connection between (in this example) 
the physics of acoustics in a fluid flow and at least some significant features 
of general relativity that we can claim to have an "analogue model of (some 
aspects of) gravity". We (and the community at large) often abuse language 
by referring to such a model as "analogue gravity" for short. 

Now the features of general relativity that one typically captures in an 
"analogue model" are the kinematic features that have to do with how fields 
(classical or quantum) are defined on curved spacetime, and the sine qua 
non of any analogue model is the existence of some "effective metric" that 
captures the notion of the curved spacetimes that arise in general relativity. 
(At the very least, one might wish to capture the notion of the Minkowski 
geometry of special relativity.) Indeed, the verbal description above (and its 
generalizations in other physical frameworks) can be converted into a precise 

lr The need for a certain degree of caution regarding the allegedly straightforward 
physics of simple fluids might be inferred from the fact that the Clay Mathematics In- 
stitute is currently offering a US$1,000,000 Millennium Prize for significant progress on 
the question of existence and uniqueness of solutions to the Navier-Stokes equation. See 
http://www.claymath.org/millennium/ for details. 

2 In correct English, the word "dumb" means "mute", as in "unable to speak". The 
word "dumb" does not mean "stupid" , though even many native English speakers get this 
wrong. 
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Figure 3: A moving fluid will drag sound pulses along with it. 

mathematical and physical statement, which ultimately is the reason that 
analogue models are of physical interest. The analogy works at two levels; 

• Geometrical acoustics. 

• Physical acoustics. 

The advantage of geometrical acoustics is that the derivation of the precise 
mathematical form of the analogy is so simple as to be almost trivial, and 
that the derivation is extremely general. The disadvantage is that in the 
geometrical acoustics limit one can deduce only the causal structure of the 
spacetime, and does not obtain a unique effective metric . The advantage of 
physical acoustics is that while the derivation of the analogy holds in a more 
restricted regime, the analogy can do more for you in that it can now specify 
a specific effective metric and accommodate a wave equation for the sound 
waves. 

2.2 Geometrical acoustics 

At the level of geometrical acoustics we need only assume that: 
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• The speed of sound c, relative to the fluid, is well defined. 

• The velocity of the fluid v, relative to the laboratory, is well defined. 

Then, relative to the laboratory, the velocity of a sound ray propagating, 
with respect to the fluid, along the direction defined by the unit vector n is 



dx 

dt 



cn + v. 



(1) 



This defines a sound cone in spacetime given by the condition n 2 



1, i.e., 



-c 2 dt 2 + (dx- vdfT 







(2) 



That is 

- [c 2 - v 2 } dt 2 - 2v • dx dt + dx ■ dx = (3) 

Solving this quadratic equation for dx as a function of dt provides a 
double cone associate with each point in space and time. This is associated 
with a conformal class of Lorentzian metrics ||378| , |388| , |390| , |391| , |286 



g' 



' -(c 2 - v 2 ) 


T " 

—V 


—v 


I 



(4) 



where Q is an unspecified but non-vanishing function. 

The virtues of the geometric approach are its extreme simplicity and 
the fact that the basic structure is dimension-independent. Moreover this 
logic rapidly (and relatively easily) generalizes to more complicated physical 



situations/ 



2.3 Physical acoustics 

It is well known that for a static homogeneous inviscid fluid the propagation 
of sound waves is governed by the simple wave equation j 22|, |225| , [266| , | 



d 1 ^ = c 2 V 2 0. (5) 



3 For instance, whenever one has a system of PDEs that can be written in first-order 
quasi-linear symmetric hyperbolic form, then it is an exact non-pcrturbative result that 
the matrix of coefficients for the first-derivative terms can be used to construct a conformal 
class of metrics that encodes the causal structure of the system of PDEs. For barotropic 
hydrodynamics this is briefly discussed in |81j| . This analysis is related to the behaviour 
of characteristics of the PDEs, and ultimately can be linked back to the Fresnel equation 
that appears in the eikonal limit. 
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Subsonic Sonic Supersonic 



1 ► 

X 

Figure 4: A moving fluid will tip the "sound cones" as it moves. Supersonic 
flow will tip the sound cones past the vertical. 

Generalizing this result to a fluid that is non- homogeneous, or to a fluid that 
is in motion, possibly even in non-steady motion, is more subtle than it at 
first would appear. To derive a wave equation in this more general situation 
we shall start by adopting a few simplifying assumptions to allow us to derive 
the following theorem. 

Theorem: If a fluid is barotropic and inviscid, and the flow is irrotational 
(though possibly time dependent) then the equation of motion for the velocity 
potential describing an acoustic disturbance is identical to the d'Alembertian 
equation of motion for a minimally coupled massless scalar field propagating 
in a (3 + l)-dimensional Lorentzian geometry 

A0 = -Lfy <T d v( j>) = 0. (6) 
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Under these conditions, the propagation of sound is governed by an acous- 
tic metric — g^ u (t,x). This acoustic metric describes a (3 + l)-dimensional 
Lorentzian (pseudo-Riemannian) geometry. The metric depends algebraically 
on the density, velocity of flow, and local speed of sound in the fluid. Specif- 
ically 



P 

c 



- C 



(7) 



(Here I is the 3x3 identity matrix.) In general, when the fluid is non- 
homogeneous and flowing, the acoustic Riemann tensor associated with this 
Lorentzian metric will be nonzero. <0 



Comment: It is quite remarkable that even though the underlying fluid dy- 
namics is Newtonian, nonrelativistic, and takes place in flat space plus time, 
the fluctuations (sound waves) are governed by a curved (3 + l)-dimensional 
Lorentzian (pseudo-Riemannian) spacetime geometry. For practitioners of 
general relativity this observation describes a very simple and concrete phys- 
ical model for certain classes of Lorentzian spacetimes, including (as we shall 
later see) black holes. On the other hand, this discussion is also potentially of 
interest to practitioners of continuum mechanics and fluid dynamics in that it 
provides a simple concrete introduction to Lorentzian differential geometric 
techniques. 



Proof: The fundamental equations of fluid dynamics | 22], [225| , |266| , 355 
are the equation of continuity 

d tP + V • (p v) = 0, (8) 

and Euler's equation (equivalent to F = m a applied to small lumps of fluid) 

P^ = P [fa + (v ■ V)v] = f . (9) 

Start the analysis by assuming the fluid to be inviscid (zero viscosity), with 
the only forces present being those due to pressure. 4 Then for the force 
density we have 

f = -Vp. (10) 

4 It is straightforward to add external forces, at least conservative body forces such as 
Newtonian gravity. 
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Via standard manipulations the Euler equation can be rewritten as 



3 t v = v X (V X v) - ~ Vp - V Q v 2 J . (11) 

Now take the flow to be vorticity free, that is, locally irrotational. Introduce 
the velocity potential such that v = — V0, at least locally. If one further 
takes the fluid to be barotropic (this means that p is a function of p only), it 
becomes possible to define 

h(p)= / -7^; so that Vh = - Vp. (12) 
Jo P(P) P 

Thus the specific enthalpy, h(p), is a function of p only. Euler's equation now 
reduces to 

-d t <P + h+^(V<P) 2 = 0. (13) 

This is a version of Bernoulli's equation. 

Now linearize these equations of motion around some assumed back- 
ground (p ,po,<j>o). Set 

P = p + epi + O(e 2 ), (14) 

P = p + e Pl + O(e 2 ), (15) 

= o + e0i + O(e 2 ). (16) 

Sound is defined to be these linearized fluctuations in the dynamical quan- 
tities. Note that this is the standard definition of [linear] sound and more 
generally of acoustical disturbances. In principle, of course, a fluid mechanic 
might really be interested in solving the complete equations of motion for the 
fluid variables (p,p, <f))- In practice, it is both traditional and extremely use- 
ful to separate the exact motion, described by the exact variables, (p,p,<p), 
into some average bulk motion, (po,Po,(fio), P ms l° w amplitude acoustic dis- 
turbances, (epi, epx, e0i). See, for example [[ ^, [225| , |266| , |355| . 



Since this is a subtle issue that we have seen cause considerable confusion 
in the past, let us be even more explicit by asking the rhetorical question: 

"How can we tell the difference between a wind gust and a sound wave?" 
The answer is that the difference is to some extent a matter of convention 

- sufficiently low-frequency long-wavelength disturbances (wind gusts) are 
conventionally lumped in with the average bulk motion. Higher-frequency, 
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shorter-wavelength disturbances are conventionally described as acoustic dis- 
turbances. If you wish to be hyper-technical, we can introduce a high-pass 
filter function to define the bulk motion by suitably averaging the exact fluid 
motion. There are no deep physical principles at stake here — merely an issue 
of convention. The place where we are making a specific physical assumption 
that restricts the validity of our analysis is in the requirement that the am- 
plitude of the high-frequency short-wavelength disturbances be small. This 
is the assumption underlying the linearization programme, and this is why 
sufficiently high-amplitude sound waves must be treated by direct solution 
of the full equations of fluid dynamics. 

Linearizing the continuity equation results in the pair of equations 

dtPo + V ■ (po v ) = 0, (17) 
dtpi + V • (pi v + po vi) = 0. (18) 

Now, the barotropic condition implies 

h{p) = h(p + ep x + 0(e 2 )) = h + e^ + 0(e 2 ). (19) 

Po 

Use this result in linearizing the Euler equation. We obtain the pair 

-«9 4 o + />o + ^(V0o) 2 = O. (20) 

-dtfa + v ± - v ■ V0i = 0. (21) 
Po 

This last equation may be rearranged to yield 

Pi = Po Wi + v - V0i). (22) 
Use the barotropic assumption to relate 

Now substitute this consequence of the linearized Euler equation into the 
linearized equation of continuity. We finally obtain, up to an overall sign, 
the wave equation: 

- d t po (<9 t 0! + v ■ V0 X ; 

Po V0! - ^ po v (d t (f)i + v ■ V0i) ) = 0. (24) 
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This wave equation describes the propagation of the linearized scalar poten- 
tial 4>i. Once 0i is determined, equation (^) determines p±, and equation 
(|23| ) then determines p\. Thus this wave equation completely determines 
the propagation of acoustic disturbances. The background fields p , p and 
v = — V0o, which appear as time- dependent and position-dependent coef- 
ficients in this wave equation, are constrained to solve the equations of fluid 
motion for a barotropic, inviscid, and irrotational flow. Apart from these 
constraints, they are otherwise permitted to have arbitrary temporal and 
spatial dependencies. 

Now, written in this form, the physical import of this wave equation is 
somewhat less than pellucid. To simplify things algebraically, observe that 
the local speed of sound is defined by 



c -2 = 

dp 

Now construct the symmetric 4x4 matrix 



(25) 



(c 2 5 ij - v l v J ) 



(26) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Then, 
introducing (3 + 1) -dimensional space-time coordinates, which we write as 
= (i; x 1 ) , the above wave equation (|24|) is easily rewritten as 



0. 



(27) 



This remarkably compact formulation is completely equivalent to equation 
(|2~4] ) and is a much more promising stepping-stone for further manipulations. 
The remaining steps are a straightforward application of the techniques of 
curved space (3 + l)-dimensional Lorentzian geometry. 

Now in any Lorentzian (that is, pseudo-Riemannian) manifold the curved 
space scalar d'Alembertian is given in terms of the metric <? MI ,(t,x) by (see, 
for example, p^, |6|, [571], |67] |167[ p3j ) 



A0 



-g 9 



(28) 



The inverse metric, g^ v {t, x), is pointwise the matrix inverse of g^ft, x), while 



g = det(g fll/ ). Thus one can rewrite the physically derived wave equation (24) 
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in terms of the d'Alembertian provided one identifies 

ST = r- (29) 

This implies, on the one hand 

detCT) = (v^) 4 g- 1 = 9- (30) 

On the other hand, from the explicit expression (|26[). expanding the deter- 
minant in minors yields 

det(r ) = (g) 4 • [(-1) • (c 2 -vl) - (-, ) 2 ] ■ [c 2 ] ■ [c 2 ] = - p l (31) 



Thus 



> = -4 



-.9 



P5 



(32) 



We can therefore pick off the coefficients of the inverse acoustic metric 



<r(*,x) = — 



1 



(33) 

We could now determine the metric itself simply by inverting this 4x4 matrix 
(and if the reader is not a general relativist, proceeding in this direct manner 
is definitely the preferred option). On the other hand, for general relativists it 
is even easier to recognize that one has in front of one a specific example of the 
Arnowitt-Deser-Misner split of a (3 + l)-dimensional Lorentzian spacetime 
metric into space + time, more commonly used in discussing initial value 
data in general relativity (See, for example, [[26711 PP 505-508.) The acoustic 
metric is then read off by inspection 



9y.u = — 



Po 

c 



~Vn 



(34) 



o • w y 

Equivalently, the acoustic interval can be expressed as 

ds 2 = dx^ dx v = ^ [-c 2 dt 2 + (dx* - vl dt) Sij {dx j - v J Q dt)} . (35) 

This completes the proof of the theorem. We have presented the theorem 
and proof, which closely follows the discussion in | || . in considerable detail 
because it is a standard template that can be readily generalised in many 
ways. This discussion can then be used as a starting point to initiate the 
analysis of numerous and diverse physical models. 
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2.4 General features of the acoustic metric 



A few brief comments should be made before proceeding further: 

• Observe that the signature of this effective metric is indeed (—,+,+,+), 
as it should be to be regarded as Lorentzian. 

• Observe that in physical acoustics it is the inverse metric density, 



that is of more fundamental significance for deriving the wave equation 
than is the metric itself. (This observation continues to hold in more 
general situations where it is often significantly easier to calculate the 
tensor density f^ u than it is to calculate the effective metric g^ u .) 

• It should be emphasized that there are two distinct metrics relevant to 
the current discussion: 

— The physical spacetime metric is in this case just the usual flat 
metric of Minkowski space: 



(Here c^t = speed of light in vacuum.) The fluid particles cou- 
ple only to the physical metric r]^ u . In fact the fluid motion is 
completely non-relativistic, so that \\v Q \\ <C ci ig ht, and it quite 
sufficient to consider Galilean relativity for the underlying fluid 
mechanics. 

— Sound waves on the other hand, do not "see" the physical metric 
at all. Acoustic perturbations couple only to the effective acoustic 
metric g^ u . 

• The geometry determined by the acoustic metric does however inherit 
some key properties from the existence of the underlying flat physical 
metric. For instance, the topology of the manifold does not depend 
on the particular metric considered. The acoustic geometry inherits 
the underlying topology of the physical metric — ordinary 3ft 4 — with 
possibly a few regions excised (due to whatever hard-wall boundary 
conditions one might wish to impose on the fluid). In systems con- 
strained to have effectively less than 3 spacelike dimensions one can 




(36) 



T]„ v = (diagf-cjjght,!,!, 



1]) 



(37) 
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reproduce more complicated topologies (consider for example an effec- 
tively one-dimensional flow in a tubular ring). 



Furthermore, the acoustic geometry automatically inherits from the un- 
derlying Newtonian time parameter, the important property of "stable 
causality" p7j p3p. Note that 



gT (V M t) (V„*) = -— < 0. (38) 
Po c 

This precludes some of the more entertaining causality-related patholo- 
gies that sometimes arise in general relativity (For a general discussion 
of causal pathologies in general relativity, see for example [|167| , |166| , |162| , 
p|, [163] g9j). 



Other concepts that translate immediately are those of "ergo-region" , 
"trapped surface", "apparent horizon", and "event horizon". These 
notions will be developed more fully in the following subsection. 

The properly normalized four-velocity of the fluid is 

V» = t^l, (39) 

so that 

9 ia ,VV = g(y,V) = -l. (40) 

This four-velocity is related to the gradient of the natural time param- 
eter by 

(\ ■ v l \ V>* 
V,t = (1,0,0,0); W = -±L-!» = --!=. (41) 

Po c VPoc 

Thus the integral curves of the fluid velocity field are orthogonal (in the 
Lorentzian metric) to the constant time surfaces. The acoustic proper 
time along the fluid flow lines (streamlines) is 



4p^~c dt, (42) 

and the integral curves are geodesies of the acoustic metric if and only 
if po c is position independent. 
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• Observe that in a completely general (3 + l)-dimensional Lorentzian 
geometry the metric has 6 degrees of freedom per point in spacetime. 
(4x4 symmetric matrix =>- 10 independent components; then subtract 
4 coordinate conditions). 

In contrast, the acoustic metric is more constrained. Being specified 
completely by the three scalars (fio(t, x), p (t, x), and c(t, x), the acous- 
tic metric has at most 3 degrees of freedom per point in spacetime. 
The equation of continuity actually reduces this to 2 degrees of free- 
dom, which can be taken to be </>o(t,x) and c(£,x). 

Thus the simple acoustic metric of this section can at best reproduce 
some subset of the generic metrics of interest in general relativity. 

• A point of notation: Where the general relativist uses the word "sta- 
tionary" the fluid dynamicist uses the phrase "steady flow". The 
general-relativistic word "static" translates to a rather messy constraint 
on the fluid flow (to be discussed more fully below). 

• Finally, we should emphasise that in Einstein gravity the spacetime 
metric is related to the distribution of matter by the non-linear Einstein- 
Hilbert differential equations. In contrast, in the present context, the 
acoustic metric is related to the distribution of matter in a simple al- 
gebraic fashion. 

2.5 Dumb holes — ergoregions, horizons, and surface 
gravity 

Let us start with the notion of an ergo-region: Consider integral curves of 
the vector 

K* = (d/dty = (1, 0, 0, 0)". (43) 

If the flow is steady then this is the time translation Killing vector. Even if 
the flow is not steady the background Minkowski metric provides us with a 
natural definition of "at rest" . Then 5 

9,u {d/dtY (d/dty = g tt = -[c 2 - v 2 }. (44) 

5 Henceforth, in the interests of notational simplicity, we shall drop the explicit subscript 
on background field quantities unless there is specific risk of confusion. 
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This quantity changes sign when | |v| | > c. Thus any region of supersonic flow 
is an ergo-region. (And the boundary of the ergo-region may be deemed to be 
the er go-surf ace.) The analogue of this behaviour in general relativity is the 
ergosphere surrounding any spinning black hole — it is a region where space 



'moves" with superluminal velocity relative to the fixed stars ||267 , |167| , |4 



A trapped surface in acoustics is defined as follows: take any closed two- 
surface. If the fluid velocity is everywhere inward-pointing and the normal 
component of the fluid velocity is everywhere greater than the local speed 
of sound, then no matter what direction a sound wave propagates, it will 
be swept inward by the fluid flow and be trapped inside the surface. The 
surface is then said to be outer-trapped. (For comparison with the usual sit- 
uation in general relativity see |167| , pages 319-323] or ||423| , pages 310-311].) 
Inner-trapped surfaces (anti-trapped surfaces) can be defined by demanding 
that the fluid flow is everywhere outward-pointing with supersonic normal 
component. It is only because of the fact that the background Minkowski 
metric provides a natural definition of "at rest" that we can adopt such a 
simple and straightforward definition. In ordinary general relativity we need 
to develop considerable additional technical machinery, such as the notion 
of the "expansion" of bundles of ingoing and outgoing null geodesies, before 
defining trapped surfaces. That the above definition for acoustic geometries 
is a specialization of the usual one can be seen from the discussion on pages 
262-263 of Hawking and Ellis ||167|| . The acoustic trapped region is now 



defined as the region containing outer trapped surfaces, and the acoustic (fu- 
ture) apparent horizon as the boundary of the trapped region. (We can also 
define anti-trapped regions and past apparent horizons but these notions are 
of limited utility in general relativity.) 6 

The event horizon (absolute horizon) is defined, as in general relativity, 
by demanding that it be the boundary of the region from which null geodesies 
(phonons) cannot escape. This is actually the future event horizon. A past 
event horizon can be defined in terms of the boundary of the region that 
cannot be reached by incoming phonons — strictly speaking this requires us 
to define notions of past and future null infinities, but we will simply take 
all relevant incantations as understood. In particular the event horizon is a 
null surface, the generators of which are null geodesies. 



6 This discussion naturally leads us to what is perhaps the central question of analogue 
models — just how much of the standard "laws of black hole mechanics" |2l], [424| carry 
over into these analogue models? Quite a lot but not everything — that's our main topic 
for the rest of the review. 
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Figure 5: A moving fluid can form "trapped surfaces" when supersonic How 
tips the sound cones past the vertical. 

In all stationary geometries the apparent and event horizons coincide, and 
the distinction is immaterial. In time-dependent geometries the distinction 
is often important. When computing the surface gravity we shall restrict 
attention to stationary geometries (steady flow). In fluid flows of high sym- 
metry, (spherical symmetry, plane symmetry) the ergosphere may coincide 
with the acoustic apparent horizon, or even the acoustic event horizon. This 
is the analogue of the result in general relativity that for static (as opposed to 
stationary) black holes the ergosphere and event horizon coincide. For many 
more details, including appropriate null coordinates and Carter-Penrose di- 
agrams, both in stationary and time-dependent situations, see ||13|| . 

Because of the definition of event horizon in terms of phonons (null 
geodesies) that cannot escape the acoustic black hole, the event horizon is 
automatically a null surface, and the generators of the event horizon are au- 
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tomatically null geodesies. In the case of acoustics there is one particular 
parameterization of these null geodesies that is "most natural" , which is the 
parameterization in terms of the Newtonian time coordinate of the under- 
lying physical metric. This allows us to unambiguously define a "surface 
gravity" even for non-stationary (time-dependent) acoustic event horizons, 
by calculating the extent to which this natural time parameter fails to be 
an affine parameter for the null generators of the horizon. (This part of the 
construction fails in general relativity where there is no universal natural 
time-coordinate unless there is a timelike Killing vector — this is why ex- 
tending the notion of surface gravity to non-stationary geometries in general 
relativity is so difficult.) 

When it comes to explicitly calculating the surface gravity in terms of 
suitable gradients of the fluid flow, it is nevertheless very useful to limit at- 
tention to situations of steady flow (so that the acoustic metric is stationary). 
This has the added bonus that for stationary geometries the notion of "acous- 
tic surface gravity" in acoustics is unambiguously equivalent to the general 
relativity definition. It is also useful to take cognizance of the fact that the 
situation simplifies considerably for static (as opposed to merely stationary) 
acoustic metrics. 

To set up the appropriate framework, write the general stationary acoustic 
metric in the form 

P 
c 



ds 2 = p - [-c 2 dt 2 + (dx - v dt) 2 } . (45) 



The time translation Killing vector is simply = (1; ), with 

K 2 = g^K»K u = -||K|| 2 = - P -\(? - v 2 }. (46) 
The metric can also be written as 

ds 2 = P r_( c 2 _ v 2) df 2 _ 2v ■ dx dt + (dx) 2 l . (47) 

c 

Now suppose that the vector v/(c 2 — v 2 ) is integrable, then we can define a 
new time coordinate by 

v ■ dx 

dr = dt + -. 48 

Cr — V z 

Substituting this back into the acoustic line element gives 



d. 2 = ^ 



c 



{ 6* + } 



(49) 
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In this coordinate system the absence of the time-space cross-terms makes 
manifest that the acoustic geometry is in fact static (the Killing vector is 
hypersurface orthogonal). The condition that an acoustic geometry be static, 
rather than merely stationary, is thus seen to be 

V x {j^} = 0, (50) 

that is, (since in deriving the existence of the effective metric we have already 
assumed the fluid to be irrotational), 

v x V(c 2 - v 2 ) = 0. (51) 

This requires the fluid flow to be parallel to another vector that is not quite 
the acceleration but is closely related to it. (Note that, because of the vortic- 
ity free assumption, | Vt> 2 is just the three-acceleration of the fluid, it is the 
occurrence of a possibly position dependent speed of sound that complicates 
the above.) 

Once we have a static geometry, we can of course directly apply all of the 
standard tricks [|374|| for calculating the surface gravity developed in general 



relativity. We set up a system of fiducial observers (FIDOS) by properly 
normalizing the time-translation Killing vector 

K K 

vfido = Tjj^jT - 7= , = = (52) 



y/{ P /c) [c 2 - V 2 } 



The four-acceleration of the FIDOS is defined as 

Afido = (V FIDO • V)V FIDO , (53) 

and using the fact that K is a Killing vector, it may be computed in the 
standard manner 

1V II K H 2 

Afido = +^ , " , ' • (54) 



That is 



1 

'-FIDO — ^ 



V(c*-v 2 ) V(p/ C ) 1 
(c 2 - v 2 ) (p/c) \ ' 1 ' 

The surface gravity is now defined by taking the norm 1 1 Afido 1 1 ; multiplying 
by the lapse function, ||K|| = \J{p/c) [c 2 — v 2 ], and taking the limit as one 
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approaches the horizon: \v\ — > c (remember that we are currently dealing 
with the static case). The net result is 



K|| = -v • V(c 2 - v 2 ) + 0(c 2 - v 2 ), (56) 



I |^-FIDO 

Z 

so that the surface gravity is given in terms of a normal derivative by 7 



= 1 d(c 2 - v 2 ) = c d(c - v) 
2 dn dn 



This is not quite Unruh's result [B7R |379|, |38(J| since he implicitly took the 



speed of sound to be a position-independent constant. The fact that prefac- 
tor p/c drops out of the final result for the surface gravity can be justified 
by appeal to the known conformal invariance of the surface gravity ||194| . 



Though derived in a totally different manner, this result is also compatible 
with the expression for "surface-gravity" obtained in the solid-state black 
holes of Reznik |[321|1 , wherein a position dependent (and singular) refractive 



index plays a role analogous to the acoustic metric. As a further consis- 
tency check, one can go to the spherically symmetric case and check that 



this reproduces the results for "dirty black holes" enunciated in [ 387 . 

Since this is a static geometry, the relationship between the Hawking 
temperature and surface gravity may be verified in the usual fast-track man- 
ner — using the Wick rotation trick to analytically continue to Euclidean 
space ||149|| . If you don't like Euclidean signature techniques (which are in 



any case only applicable to equilibrium situations) you should go back to the 
original Hawking derivations [164, |165|| . 8 



One final comment to wrap up this section: the coordinate transform we 
used to put the acoustic metric into the explicitly static form is perfectly good 
mathematics, and from the general relativity point of view is even a simplifi- 
cation. However, from the point of view of the underlying Newtonian physics 
of the fluid, this is a rather bizarre way of deliberately de-synchronizing your 
clocks to take a perfectly reasonable region — the boundary of the region of 
supersonic flow — and push it out to "time" plus infinity. From the fluid 
dynamics point of view this coordinate transformation is correct but per- 
verse, and it is easier to keep a good grasp on the physics by staying with 
the original Newtonian time coordinate. 

7 Because of the background Minkowski metric there can be no possible confusion as to 
the definition of this normal derivative. 

8 There are a few potential subtleties in the derivation of the existence Hawking radia- 
tion which we are for the time being glossing over, see Section for details. 



26 



If the fluid flow does not satisfy the integrability condition which allows us 
to introduce an explicitly static coordinate system, then defining the surface 
gravity is a little trickier. 

Recall that by construction the acoustic apparent horizon is in general 
defined to be a two-surface for which the normal component of the fluid 
velocity is everywhere equal to the local speed of sound, whereas the acoustic 
event horizon (absolute horizon) is characterized by the boundary of those 
null geodesies (phonons) that do not escape to infinity. In the stationary 
case these notions coincide, and it is still true that the horizon is a null 
surface, and that the horizon can be ruled by an appropriate set of null curves. 
Suppose we have somehow isolated the location of the acoustic horizon, then 
in the vicinity of the horizon we can split up the fluid flow into normal and 
tangential components 

v = v_l + vn; where vj_ = v± n. (58) 

Here (and for the rest of this particular section) it is essential that we use 
the natural Newtonian time coordinate inherited from the background New- 
tonian physics of the fluid. In addition n is a unit vector field that at the 
horizon is perpendicular to it, and away from the horizon is some suitable 
smooth extension. (For example, take the geodesic distance to the horizon 
and consider its gradient.) We only need this decomposition to hold in some 
open set encompassing the horizon and do not need to have a global de- 
composition of this type available. Furthermore, by definition we know that 
vj_ = c at the horizon. Now consider the vector field 

L" = (l; vjj). (59) 

Since the spatial components of this vector field are by definition tangent to 
the horizon, the integral curves of this vector field will be generators for the 
horizon. Furthermore the norm of this vector (in the acoustic metric) is 



(c 2 — v 2 ) — 2v|| • v + vy ■ V|| 



P -{c 2 -vl). (60) 

c 



In particular, on the acoustic horizon L M defines a null vector field, the in- 
tegral curves of which are generators for the acoustic horizon. We shall now 
verify that these generators are geodesies, though the vector field L is not 
normalized with an afline parameter, and in this way shall calculate the 
surface gravity. 
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Consider the quantity (L • V)L and calculate 



To calculate the first term note that 

P < r 2 2 



Thus 



-v, [f i 



c 2 — -u 2 



(f ^) [H 



And so: 



vii ■ V 



P I 2 2 
LC (C 



P / 2 2 
-(C -T, ± 



J] 



( - i; x 



On the horizon, where c = v±, this simplifies tremendously 

P fn.r, ,,2 P d(c 2 -vl) 



dn 



(L a L [f3M )\ hovizon = (0; V^c 2 - v\)) 
Similarly, for the second term we have 

V /3 (L 2 ) = (0;V, 
On the horizon this again simplifies 

V,(^ 2 )| h on Z on = + P ~ (0; V,(c 2 - vl)) = + P - ^ ~ ^ (0; n,) . 
c c an 

There is partial cancellation between the two terms, and so 

1 p d(c 2 — v 2 ' 



2 C 



<9n 



(0;%). 



while 



(^- / /i)horizon (0, C Tij) 



(61) 



(62) 



(63) 



• (64) 



(0;n,). (65) 



(66) 



(67) 

(68) 
(69) 
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Comparing this with the standard definition of surface gravity [[ 



\L V tt Iji) horizon (-^/i) horizon' 



(70) 



we finally have 



9h = 



1<9(c 2 
2 



) d(c-v ± ) 



(71) 



<9n <9?7. 

This is in agreement with the previous calculation for static acoustic black 
holes, and insofar as there is overlap, is also consistent with results of Un- 
ruh |78[ |79[ M , Reznik [l321j , and the results for "dirty black holes" 
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From the construction it is clear that the surface gravity is a measure of the 
extent to which the Newtonian time parameter inherited from the underlying 
fluid dynamics fails to be an affine parameter for the null geodesies on the 
horizon. 10 

Again, the justification for going into so much detail on this specific model 
is that this style of argument can be viewed as a template — it will (with 
suitable modifications) easily generalize to more complicated analogue mod- 
els. 



2.5.1 Example: Vortex geometry 

As an example of a fluid flow where the distinction between ergosphere and 
acoustic event horizon is critical consider the "draining bathtub" fluid flow. 
We shall model a draining bathtub by a (3 + 1) dimensional flow with a 
linear sink along the z axis. Let us start with the simplifying assumption 
that the background density p is a position-independent constant throughout 
the flow (which automatically implies that the background pressure p and 
speed of sound c are also constant throughout the fluid flow). The equation 
of continuity then implies that for the radial component of the fluid velocity 
we must have 



9 Thcrc is an issue of normalization here. On the one hand we want to be as close as 
possible to general relativistic conventions. On the other hand, we would like the surface 
gravity to really have the dimensions of an acceleration. The convention adopted here, 
with one explicit factor of c, is the best compromise we have come up with. (Note that 
in an acoustic setting, where the speed of sound is not necessarily a constant, we cannot 
simply set c — » 1 by a choice of units.) 

10 There are situations in which this surface gravity is a lot larger than one might naively 



expect |]241 ] 
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v f oc i. (72) 

In the tangential direction, the requirement that the flow be vorticity free 
(apart from a possible delta-function contribution at the vortex core) implies, 
via Stokes' theorem, that 

u*oc~ (73) 
r 

(If these flow velocities are nonzero, then following the discussion of [ |102| 



there must be some external force present to set up and maintain the back- 
ground flow. Fortunately it is easy to see that this external force affects only 
the background flow and does not influence the linearized fluctuations we are 
interested in.) 

For the background velocity potential we must then have 

0(r, 6) = -A \n(r/a) - B 9. (74) 

Note that, as we have previously hinted, the velocity potential is not a true 
function (because it has a discontinuity on going through 2ir radians). The 
velocity potential must be interpreted as being defined patch-wise on over- 
lapping regions surrounding the vortex core at r = 0. The velocity of the 
fluid flow is 

„, (Ar + B§) 
v=-V0=- '-. 75 

r 

Dropping a position-independent prefactor, the acoustic metric for a 
draining bathtub is explicitly given by 

ds 2 = -c 2 dt 2 + (dr - jd?j + (rd9- ^-d?j + dz 2 . (76) 
Equivalent ly 

ds 2 = - (c 2 - A2 ^ 2 B2 ^J dt 2 -2j dr dt-2B d6 dt + dr 2 + r 2 d6 2 + dz 2 . (77) 

A similar metric, restricted to A = (no radial flow), and generalized to an 
anisotropic speed of sound, has been exhibited by Volovik ||405|| , that metric 
being a model for the acoustic geometry surrounding physical vortices in 
superfluid 3 He. (For a survey of the many analogies and similarities between 
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Figure 6: A collapsing vortex geometry (draining bathtub): The green spi- 
rals denote streamlines of the fluid flow. The outer circle represents the 
ergosurface while the inner circle represents the [outer] event horizon. 



the physics of superfluid 3 He and the Standard Electroweak Model see ||421|| , 
this reference is also useful as background to understanding the Lorentzian 
geometric aspects of 3 He fluid flow.) Note that the metric given above is not 
identical to the metric of a spinning cosmic string, which would instead take 



the form 11389 



ds 2 = -c 2 (di - Add) 2 + dr 2 + (1 - B)r 2 dd 2 + dz 2 . 



(78) 



In conformity with previous comments, the vortex fluid flow is seen to possess 
an acoustic metric that is stably causal and which does not involve closed 
timelike curves. (At large distances it is possible to approximate the vortex 
geometry by a spinning cosmic string || . | . but this approximation becomes 
progressively worse as the core is approached.) 
The ergosphere forms at 



vM 2 + B 2 



ergosphere 



(79) 
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Note that the sign of A is irrelevant in defining the ergosphere and ergo- 
region: it does not matter if the vortex core is a source or a sink. 

The acoustic event horizon forms once the radial component of the fluid 
velocity exceeds the speed of sound, that is at 

'"horizon • (^0) 

The sign of A now makes a difference. For A < we are dealing with a 
future acoustic horizon (acoustic black hole), while for A > we are dealing 
with a past event horizon (acoustic white hole). 



2.5.2 Example: Slab geometry 

A popular model for the investigation of event horizons in the acoustic anal- 
ogy is the one-dimensional slab geometry where the velocity is always along 
the z direction and the velocity profile depends only on z. The continuity 
equation then implies that p(z) v(z) is a constant, and the acoustic metric 
becomes 

ds 2 oc , r~ . r \-c{z) 2 dt 2 + {dz - v{z) dt} 2 + dx 2 + dy 2 } . (81) 
v{z) c(z) 

That is 

ds 2 oc , } , . [- \c(z) 2 - viz) 2 } dt 2 - 2v(z) dz dt + dx 2 + dy 2 + dz 2 ] . 
v(z) C[Z) 

(82) 

If we set c = 1 and ignore the conformal factor we have the toy model 
acoustic geometry discussed by Unruh ||38C| , page 2828, equation (8)] Jacob- 
190| , page 7085, equation (4)], Corley and Jacobson fl90"[ , and Corley f 



son 



(In this situation one must again invoke an external force to set up and main- 
tain the fluid flow. Since the conformal factor is regular at the event horizon, 
we know that the surface gravity and Hawking temperature are independent 



of this conformal factor ||194j| .) In the general case it is important to realise 
that the flow can go supersonic for either of two reasons: the fluid could speed 
up, or the speed of sound could decrease. When it comes to calculating the 
"surface gravity" both of these effects will have to be taken into account. 
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2.5.3 Example: Painleve-Gullstrand geometry 

To see how close the acoustic metric can get to reproducing the Schwarzschild 
geometry it is first useful to introduce one of the more exotic representations 
of the Schwarzschild geometry: the Painleve-Gullstrand line element, which 
is simply an unusual choice of coordinates on the Schwarzschild spacetime. 11 
In modern notation the Schwarzschild geometry in ingoing (+) and outgoing 
(— ) Painleve-Gullstrand coordinates may be written as: 



ds 2 = -dt 2 



dr ± 



2GM 



-dt 



(d6 2 + sin 2 9 d<j) 2 ) . 



(83) 



Equivalently 



ds' 



1 



r I 



2GM 



drdt + dr 2 + r 2 (d8 2 + sin 2 6 d0 2 ) . (84) 



This representation of the Schwarzschild geometry was not (until the 
advent of the analogue models) particularly well-known, and it has been 
independently rediscovered several times during the 20th century. See for in- 
stance Painleve ||295|| , Gullstrand ||156|| , Lemaitre ||230|| , the related discussion 
by Israel [|168jl , and more recently, the paper by Kraus and Wilczek |221| . The 
Painleve-Gullstrand coordinates are related to the more usual Schwarzschild 
coordinates by 



tpc = ts± 



4M arctanh 



'2GM 



2 V2GMr 



Or equivalently 



dt 



PG 



dt,± 



dr. 



(85) 



(86) 



1 - 2GM/r 

With these explicit forms in hand, it becomes an easy exercise to check the 
equivalence between the Painleve-Gullstrand line element and the more usual 
Schwarzschild form of the line element. It should be noted that the + sign 
corresponds to a coordinate patch that covers the usual asymptotic region 
plus the region containing the future singularity of the maximally extended 
Schwarzschild spacetime. It thus covers the future horizon and the black hole 



11 The Painleve-Gullstrand line clement is sometimes called the Lemaitre line clement. 
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singularity. On the other hand the — sign corresponds to a coordinate patch 
that covers the usual asymptotic region plus the region containing the past 
singularity. It thus covers the past horizon and the white hole singularity. 

As emphasized by Kraus and Wilczek, the Painleve-Gullstrand line el- 
ement exhibits a number of features of pedagogical interest. In particular 
the constant time spatial slices are completely flat — the curvature of space 
is zero, and all the spacetime curvature of the Schwarzschild geometry has 
been pushed into the time-time and time-space components of the metric. 

Given the Painleve-Gullstrand line element, it might seem trivial to force 
the acoustic metric into this form: simply take p and c to be constants, and 
set v = ^2GM/rl While this certainly forces the acoustic metric into the 
Painleve-Gullstrand form the problem with this is that this assignment is 
incompatible with the continuity equation V ■ (pv) 7^ that was used in 
deriving the acoustic equations. 

The best we can actually do is this: Pick the speed of sound c to be a 
position independent constant, which we normalize to unity (c = 1). Now 
set v = \j2GMjr, and use the continuity equation V • (pv) = to deduce 
p|v| oc 1/r 2 so that p oc r~ 3 / 2 . Since the speed of sound is taken to be 
constant we can integrate the relation c 2 = dp/ dp to deduce the equation 
of state must be p = p^ + c 2 p and that the background pressure satisfies 
p — Poo oc c 2 r~ 3 / 2 . Overall the acoustic metric is now 



ds 2 oc r" 3 / 2 



-dt 2 +[dr±J 2 —dt ) + r 2 (dd 2 + sin 2 9 d<f) 2 ) 



(87) 



So we see that the net result is conformal to the Painleve-Gullstrand 
form of the Schwarzschild geometry but not identical to it. For many pur- 
poses this is quite good enough: we have an event horizon, we can define 
surface gravity, we can analyze Hawking radiation. 12 Since surface gravity 



and Hawking temperature are conformal invariants ||194|| this is sufficient for 
analyzing basic features of the Hawking radiation process. The only way in 
which the conformal factor can influence the Hawking radiation is through 
backscattering off the acoustic metric. (The phonons are minimally coupled 
scalars, not conformally coupled scalars so there will in general be effects on 
the frequency-dependent greybody factors.) 



12 Similar constructions work for the Reissner-Nordstrom geometry [241[, as long as one 
does not get too close to the sing ularity. Likewise certain aspects of the Kerr geometry 
can be emulated in this way [|402| . 
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If we focus attention on the region near the event horizon, the conformal 
factor can simply be taken to be a constant, and we can ignore all these 
complications. 



2.6 Regaining geometric acoustics 

Up to now, we have been developing general machinery to force acoustics into 
Lorentzian form. This can be justified either with a view to using fluid me- 
chanics to teach us more about general relativity, or to using the techniques 
of Lorentzian geometry to teach us more about fluid mechanics. 

For example, given the machinery developed so far, taking the short 
wavelength/ high frequency limit to obtain geometrical acoustics is now easy. 
Sound rays (phonons) follow the null geodesies of the acoustic metric. Com- 
pare this to general relativity where in the geometrical optics approximation 
light rays (photons) follow null geodesies of the physical spacetime metric. 
Since null geodesies are insensitive to any overall conformal factor in the met- 
ric ||2671 , |167| , [423|| one might as well simplify life by considering a modified 



conformally related metric 



-{c 2 - vl) \ -v 3 



()'■' 



(88) 



This immediately implies that, in the geometric acoustics limit, sound prop- 
agation is insensitive to the density of the fluid. In this limit, acoustic prop- 
agation depends only on the local speed of sound and the velocity of the 
fluid. It is only for specifically wave related properties that the density of 
the medium becomes important. 

We can rephrase this in a language more familiar to the acoustics commu- 
nity by invoking the Eikonal approximation. Express the linearized velocity 
potential, 0i, in terms of an amplitude, a, and phase, ip, by 0i ~ ae llf . Then, 
neglecting variations in the amplitude a, the wave equation reduces to the 
Eikonal equation 

hT d u ip = 0. (89) 

This Eikonal equation is blatantly insensitive to any overall multiplicative 
prefactor (conformal factor). 

As a sanity check on the formalism, it is useful to re-derive some standard 
results. For example, let the null geodesic be parameterized by x^{t) = 
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(i;x(i)). Then the null condition implies 



, dx^ 1 dx v 
h^u— — = 



dt dt 



- c" 



2vl 



„■ dx l dx % dx l 
dt dt dt 



<ix 



c. 



(90) 



Here the norm is taken in the flat physical metric. This has the obvious 
interpretation that the ray travels at the speed of sound, c, relative to the 
moving medium. 

Furthermore, if the geometry is stationary one can do slightly better. Let 
x^(s) = (t(s); x(s)) be some null path from xj to x 2 , parameterized in terms 



of physical arc length (i.e. ||dx/ds|| = 1). 
path is 

dx^ ( dt dx' 
ds yds'' ds 



Then the tangent vector to the 



(91) 



The condition for the path to be null (though not yet necessarily a null 
geodesic) is 

dx^ dx v „ 
ds ds 

Using the explicit algebraic form for the metric, this can be expanded to 
show 



- c -v n 



dt 
ds 



dt 
ds 



0. 



Solving this quadratic 



(93) 



dt\ 
ds~) 



dx 1 
I ds 



C Z — Vn 



(94) 



Therefore, the total time taken to traverse the path is 



(dt/ds) ds 



Xl 



/ \]{c 2 ~ rf)ds 2 + (vl dx*) 2 - vl dx'}. (95) 
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If we now recall that extremizing the total time taken is Fermat's principle 
for sound rays, we see that we have checked the formalism for stationary 
geometries (steady flow) by reproducing the discussion on page 262 of Landau 
and Lifshitz p25| . 



2.7 Generalizing the physical model 

There are a large number of ways in which the present particularly sim- 
ple analogue model can be generalized. Obvious issues within the current 
physical framework are: 



• Adding external forces. 

• Working in (1+1) or (2+1) dimensions. 

• Adding vorticity, to go beyond the irrotational constraint. 



Beyond these immediate questions, we could also seek similar effects in other 
physical or mathematical frameworks. 



2.7.1 External forces 



Adding external forces is easy, an early discussion can be found in ||391|| and 
more details are available in [ | ■ .-|| . The key point is that with an external 
force one can to some extent shape the background flow (see for example 
the discussion on 



Upon linearization, the fluctuations are however 



insensitive to any external force. 



2.7.2 The role of dimension 

The role of spacetime dimension in these acoustic geometries is sometimes 
a bit surprising and potentially confusing. This is important because there 
is a real physical distinction, for instance, between truly (2+l)-dimensional 
systems and effectively (2+l)-dimensional systems in the form of (3+1)- 
dimensional systems with cylindrical symmetry. Similarly there is a real 
physical distinction between a truly (l+l)-dimensional system and a (3+1)- 
dimensional system with transverse symmetry. We emphasise that in carte- 
sian coordinates the wave equation 
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where 



-Pic 2 



f 



flu 



d_ 

dx h 



-p v j /c 2 



-p v l /c 2 p {5 lj — v l v j /c 2 } 



(96) 
(97) 



holds independent of the dimensionality of spacetime. It depends only on 
the Euler equation, the continuity equation, a barotropic equation of state, 



and the assumption of irrotational flow | [8|, |388| , |390| , |391|| . 
Introducing the inverse acoustic metric g^, defined by 
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the wave equation (^) corresponds to the d'Alembertian wave equation in 
a curved space-time with contravariant metric tensor: 



flV 



2/(d-l) 


" -1/c 2 


-v T /c 2 




-v/c 2 


Idxd- v®v T /c 2 



(99) 



where d is the dimension of space (not spacetime). The covariant acoustic 
metric is then 
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d = 3: The acoustic line element for three space and one time dimension 
reads 

. c, 
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d = 2: The acoustic line element for two space and one time dimension 
reads 
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—v 
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(102) 



This situation would be appropriate, for instance, when dealing with surface 
waves or excitations confined to a particular substrate. 
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d = 1: The naive form of the acoustic metric in (1 + 1) dimensions is ill- 
defined, because the conformal factor is raised to a formally infinite power. 
This is a side effect of the well-known conformal invariance of the Laplacian 
in 2 dimensions. The wave equation in terms of the densitized inverse metric 
j/j-f continues to make good sense; it is only the step from f^ u to the effective 
metric that breaks down. 

Note that this issue only presents a difficulty for physical systems that are 
intrinsically one-dimensional. A three-dimensional system with plane sym- 
metry, or a two-dimensional system with line symmetry, provides a perfectly 
well behaved model for (1+1) dimensions, as in the cases d = 3 and d = 2 
above. 



2.7.3 Adding vorticity 

For the preceding analysis to hold it is necessary and sufficient that the flow 
locally be vorticity free, V x v = 0, so that velocity potentials exist on 
an atlas of open patches. Note that the irrotational condition is automat- 
ically satisfied for the super-fluid component of physical superfluids. (This 
point has been emphasised by Comer [~]. who has also pointed out that in 



superfluids there will be multiple acoustic metrics — and multiple acoustic 
horizons — corresponding to first and second sound). Even for normal fluids, 
vorticity free flows are common, especially in situations of high symmetry. 
Furthermore, the previous condition enables us to handle vortex filaments, 
where the vorticity is concentrated into a thin vortex core, provided we do 
not attempt to probe the vortex core itself. It is not necessary for the velocity 
potential <fi to be globally defined. 

Though physically important, dealing with situations of distributed vor- 
ticity is much more difficult, and the relevant wave equation is more compli- 
cated in that the velocity scalar is now insufficient to completely characterize 
the fluid flow. 13 An approach similar to the spirit of the present discussion, 
but in terms of Clebsch potentials, can be found in [|309|| . The eikonal approx- 
imation (geometrical acoustics) leads to the same conformal class of metrics 
previously discussed, but in the realm of physical acoustics the wave equation 
is considerably more complicated than a simple d'Alembertian. (Roughly 
speaking, the vorticity becomes a source for the d'Alembertian, while the 

13 Vorticity is automatically generated, for instance, whenever the background fluid is 
non-barotropic, and in particular when Vp x Vp 7^ 0. Furthermore, it has been argued 



in [345] that quantum back-reaction can also act as a source for vorticity. 
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vorticity evolves in response to gradients in a generalized scalar potential. 
This seems to take us outside the realm of models of direct interest to the 
general relativity community.) 14 



2.8 Simple Lagrangian meta- model 

As a first (and rather broad) example of the very abstract ways in which 
the notion of an acoustic metric can be generalized, we start from the simple 
observation that irrotational barotropic fluid mechanics can be described by 
a Lagrangian, and ask if we can extend the notion of an acoustic metric to 
all (or at least some wide class of) Lagrangian systems? 

Indeed, suppose we have a single scalar field whose dynamics is governed 
by some generic Lagrangian £(5^0,0), which is some arbitrary function of 
the field and its first derivatives (here we will follow the notation and ideas 
of [0). In the general analysis that follows the previous irrotational and 
inviscid fluid system is included as a particular case; the dynamics of the 
scalar field is now much more general. We want to consider linearized 
fluctuations around some background solution o (£,x) of the equations of 
motion, and to this end we write 



0(£, x) = (f> (t, x) + e0! (t, x) + -0 2 (t, x) + 0(e 3 ). 



(103) 



Now use this to expand the Lagrangian around the classical solution 0o(£, x): 
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14 



In references 1 143 



145 



141, h 42 L 1471 146, 140 



the author has attempted to argue 
that vorticity can be related to the concept of torsion in a general afhne connexion. We 
disagree. Although deriving a wave equation in the presence of vorticity very definitely 
moves one beyond the realm of a simple Riemannian spacetime, adding torsion to the 
connexion is not sufficient to capture the relevant physics. 
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It is particularly useful to consider the action 



S[4>] 



d d+l x £(<9 



(105) 



since doing so allows us to integrate by parts. (Note that the Lagrangian 
L is taken to be a scalar density, not a true scalar.) We can now use the 
Euler-Lagrange equations for the background field 



dC 



to discard the linear terms (remember we are linearizing around a solution 
of the equations of motion) and so we get 



S[4>] = S[4> ] 



id+l 
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Having set things up this way, the equation of motion for the linearized 
fluctuation is now easily read off as 
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This is a second-order differential equation with position-dependent coeffi- 
cients (these coefficients all being implicit functions of the background field 

00). 

This can be given a nice clean geometrical interpretation in terms of a 
d'Alembertian wave equation — provided we define the effective spacetime 
metric by 

d 2 C 
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(109) 



4>vi 



Note that this is another example of a situation in which calculating the 
inverse metric density is easier than calculating the metric itself. 

Suppressing the 0o except when necessary for clarity, this implies [in 
(d+1) dimensions, d space dimensions plus 1 time dimension] 

/_ a \(d-i)/2 _ _ d t { 
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Therefore 
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And, taking the inverse 
SvOo) = ( -det 
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We can now write the equation of motion for the linearized fluctuations in 
the geometrical form 



[A(s(0o))-Wo)]0i = O, 



(113) 



where A is the d'Alembertian operator associated with the effective met- 
ric g((po), and V(4> ) is the background-field-dependent (and so in general 
position-dependent) "mass term" : 
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d(d^) d{d v <P) 
d 2 C „ f d 2 C 



115) 



Thus V(<fro) is a true scalar (not a density). Note that the differential equa- 
tion (|113|) is automatically formally self-adjoint (with respect to the measure 
•>/— p d d+1 x). 

It is important to realise just how general the result is (and where the 
limitations are): it works for any Lagrangian depending only on a single 
scalar field and its first derivatives. The linearized PDE will be hyperbolic 
(and so the linearized equations will have wave-like solutions) if and only if 
the effective metric g^ u has Lorentzian signature ±[— Observe that if 
the Lagrangian contains nontrivial second derivatives you should not be too 
surprised to see terms beyond the d'Alembertian showing up in the linearized 
equations of motion. 
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As a specific example of the appearance of effective metrics due to La- 
grangian dynamics we reiterate the fact that inviscid irrotational barotropic 
hydrodynamics naturally falls into this scheme (which is why, with hindsight, 
the derivation of the acoustic metric presented earlier in this review was so 
relatively straightforward). In inviscid irrotational barotropic hydrodynam- 
ics the lack of viscosity (dissipation) guarantees the existence of a Lagrangian; 
which a priori could depend on several fields. Since the flow is irrotational 
v = — V0 is a function only of the velocity potential, and the Lagrangian 
is a function only of this potential and the density. Finally the equation of 
state can be used to eliminate the density leading to a Lagrangian that is a 
function only of the single field <fi and its derivatives. [| 

2.9 Going further 

The class of analogue models based on fluid mechanics is now quite large 
and the literature is extensive. Most of the relevant discussion will be de- 
ferred until subsequent sections, so for the time being we shall just mention 
reasonably immediate generalizations such as: 

• Working with specific fluids. 

— Superfluids. 

— Bose-Einstein condensates. 

• Abstract generalizations. 

— Normal modes in generic systems. 

— Multiple signal speeds. 

We next turn to a brief historical discussion, seeking to place the work of the 
last decade into its proper historical perspective. 
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3 History and motivation 



From the point of view of the general relativity community the history of 
analogue models can reasonably neatly (but superficially) be divided into a 
"historical" period (essentially pre-1981) and a "modern" period (essentially 
post-1981). 



3.1 Modern period 
3.1.1 1981-1999: 

The key event in the "modern" period (though largely unrecognized at the 
time) was the 1981 publication of Unruh's paper "Experimental black hole 
evaporation" ||378|| , which implemented an analogue model based on fluid 



flow, and then used the power of that analogy to probe fundamental issues 
regarding Hawking radiation from "real" general relativity black holes. 

We believe that Unruh's 1981 article represents the first observation of 
the now widely established fact that Hawking radiation has nothing to do 
with general relativity per se, but that Hawking radiation is instead a funda- 
mental curved-space quantum field theory phenomenon that occurs whenever 
a horizon is present in an effective geometry. 15 Though Unruh's 1981 paper 
was seminal in this regard, it lay largely unnoticed for many years. 

Some 10 years later Jacobson's article "Black-hole evaporation and ul- 
trashort distances" [ |188|| used Unruh's analogy to build a physical model for 
the "trans-Planckian modes" believed to be relevant to the Hawking radi- 
ation process. Progress then sped up with the relatively rapid appearance 
of ||189|| and ||379| , [380|| . (This period also saw the independent rediscovery of 
the fluid analogue model by one of the present authors | |388| , and the first 
explicit consideration of superfluids in this regard |]86|.) 



The later 1990 's then saw continued work by Jacobson and his group i 



190| , |90| , |9~^ , [200|| , with new and rather different contributions coming in the 
form of the solid state models considered by Reznik [|321|, |320| . 16 This pe- 



riod also saw the introduction of the more general class of superfluid models 



15 We emphasise: To get Hawking radiation you need an effective geometry, a horizon, 
and a suitable quantum field theory on that geometry. 



16 Rcfcrence [176| is an attempt at connecting Hawking evaporation with the physics of 
collapsing bubbles. This was part of a more general programme aimed at connecting black 



hole thermodynamics with perfect fluid thermodynamics [177]. 
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considered by Volovik and his collaborators 



TO p|, cm m ra TO 



more precise formulations of the notions of horizon, er- 



gosphere, and surface gravity in analogue models [|39 1| , |390|| , and discussions 



of the implications of analogue models regarding Bekenstein-Hawking en- 
tropy [|392| , [39C|1 . Finally, analogue spacetimes based on special relativistic 
acoustics were considered in [ 3^ . 

By the year 2000, articles on one or another aspect of analogue gravity 
were appearing at the rate of over 20 per year, and it becomes impractical 
to summarize more than a selection of them. 



3.1.2 2000: 



Key developments in 2000 were the introduction, by Garay and collaborators 
of the use of Bose-Einstein condensates as a working fluid 





138, 139 



and 

Further afield, the 



the extension of those ideas by the present authors 
trans-Planckian problem also reared its head in the context of cosmological 
inflation, and analogue model ideas previously applied to Hawking radiation 
were reused in that context ||208| , |2 



That year also marked the appearance of a review article on superfluid 
more work on "near-horizon" physics [ 5 



analogues [■ 



and the trans- 
ference of the idea of analogue-inspired "multiple metric" theories into cos- 
mology where they can be used as the basis for a precise definition of what 
is meant by a VSL 



'variable speed of light" 



on nonlinear electrodynamics were investigated in [! 
were reconsidered in [|195| , ^ 



cosmology p8[ . Models based 
3 He-A based models 



and "slow light" models in quantum di- 



electrics were considered in [|237| , [238| , |233 



The most radical proposal to appear in 2000 was that of Laughlin et 
al. I 7j. Based on taking a superfluid analogy rather literally they mooted 



an actual physical breakdown of general relativity at the horizon of a black 
hole |7|. 

Additionally, the workshop on "Analogue models of general relativity", 
held at CBPF (Rio de Janeiro) gathered some 20 international participants 
and greatly stimulated the field, leading ultimately to the publication of the 
book |86| in 2002. 
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3.1.3 2001: 



This year saw more applications of analogue-inspired ideas to cosmological 
265] , i 



inflation f08 



210, [27711, to neutron star cores 66 , and to the 



cosmological constant J 15], f416j . 

Closer to the heart of the analogue programme were the development 
of a "normal mode" analysis in |L5], [16], |399|| , the development of dielectric 



analogues in ||344j| , speculations regarding the possibly emergent nature of 
Einstein gravity [F2(], |399| , and further developments regarding the use of 
3 He- A 



107|| as an analogue for electromagnetism. Experimental proposals 

m em 



were considered in 

Vorticity was discussed in [ |.. : j| . and the use of BECs as a model for the 
breakdown of Lorentz invariance in 



39^ 



linear electrodynamics were discussed in 
vortex were investigated in [ |122| . 



Analogue models based on non- 
10^]. Acoustics in an irrotational 



The excitation spectrum in superfluids, specifically the fermion zero modes, 
were investigated in ||413| , 



while the relationship between rotational fric- 
tion in superfluids and super-radiance in rotating spacetimes was discussed 
|57| . More work on "slow light" appeared in |4l|. The possible role of 



in 



Lorentz violations at ultra-high energy was emphasized in | 



3.1.4 2002: 



"What did we learn from studying acoustic black holes?" was the title and 
theme of Parentani's article in 2002 [|302jl , while Schutzhold and Unruh devel- 
oped a rather different fluid-based analogy based on gravity waves in shallow 
water ||346| , |347|| . Super-radiance was investigated in 
gation of phonons and quasiparticles was discussed in 
on "slow light" appeared in ||126| , [313|| . 

The stability of an acoustic white hole was investigated in |[236||, while 



27 



while the propa- 
12311. More work 



further developments regarding analogue models based on nonlinear electro- 
dynamics were presented by Novello and collaborators in [|103j , |104|, |284|, [280| , 
|128||. Analogue spacetimes relevant to braneworld cosmologies were consid- 



ered in 12 



Though analogue models lead naturally to the idea of high-energy vio- 
lations of Lorentz invariance, it must be stressed that definite observational 
evidence for violations of Lorenz invariance is lacking — in fact there are 
rather strong constraints on how strong any possible Lorentz violating effect 
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might be JT37], [Hj6 . 



3.1.5 2003: 



That year saw further discussion of analogue-inspired models for black hole 
entropy and the cosmological constant ||420| , |422|| , and the development of 



116, 115, 17, 106 



There were 



analogue models for FRW geometries 
several further developments regarding the foundations of BEC-based mod- 
els in |18| , |117|| , while analogue spacetimes in superfluid neutron stars were 
further investigated in [|67 |. 

Effective geometry was the theme in 



while applications of nonlin- 
ear electrodynamics (and its effective metric) to cosmology were presented 
83f . Super-radiance was further investigated in [26], 



m 



while the limita- 
tions of the "slow light" analogue were explained in ||381| . Vachaspati argued 
for an analogy between phase boundaries and acoustic horizons in ||383 . 
Emergent relativity was again addressed in [[229| . 

The review article by Burgess |53| , emphasised the role of general rel- 
ativity as an effective field theory — the sine qua non for any attempt at 
interpreting general relativity as an emergent theory. The lecture notes by 



Jacobson [ 193 1 give a nice introduction to Hawking radiation and its connec- 
tion to analogue spacetimes. 



3.1.6 2004: 

The year 2004 saw the appearance of some 30 articles on (or closely related 
to) analogue models. Effective geometries in astrophysics were discussed by 
Perez-Bergliaffa [[30811 , while the physical realizability of acoustic Hawking 

|3g4[ 



radiation was addressed in 



in B384 425 



was invoked in 98 



More cosmological issues were raised 
while a specifically astrophysical use of the acoustic analogy 

SB 



BEC-based horizons were again considered in ||151| , |1 50|1 , while backreac- 
tion effects were the focus of attention in [JR], ^ BTE . More issues relating 
to the simulation of FRW cosmologies were raised in ||1 20| , 



121 



Unruh and Schutzhold discussed the universality of the Hawking effect PS1 



and a new proposal for possibly detecting Hawking radiation in a electro- 
magnetic wave guide [|349|1 . The causal structure of analogue spacetimes was 
considered in | I3| , while quasinormal modes attracted attention in [[31 
|&H p71|l - Two dimensional analogue models were considered in | 



239 
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There were attempts at modelling the Kerr geometry [ U)2 |. a proposal for 
using analogue models to generate massive phonon modes in BECs ||401|| , and 
an extension of the usual formalism for representing weak-field gravitational 



lensing in terms of an analogue refractive index [38 



Finally we mention the development of yet more strong observational 
bounds on possible ultra high energy Lorentz violation 



198, 199 



3.1.7 2005: 

The first few months of 2005 have seen continued and vigourous activity on 
the analogue model front. 

More studies of the super-resonance phenomenon have appeared ]25j, |114j , 
[212| , |356|| , and a mini-survey was presented in |J3J . Quasinormal modes have 
again received attention in 



79 



while the Magnus force is reanalyzed in terms 
of the acoustic geometry in \' 3|. Singularities in the acoustic geometry are 
considered in 



56 , while back-reaction has received more attention in 345 



Interest in analogue models is intense and shows no signs of abating. 



3.1.8 Analogue spacetimes as background gelstat 

To further complicate matters, there is large body of work for which analogue 
spacetime ideas provide part of the background gelstat, even if the specific 
connection may sometimes be somewhat tenuous. Among such articles we 
mention |, g | |, |7|, ||, |3] ||, |0|, [|| |3j, |37|, [|| [6(| |l], H, [68|, gjj | 



3 



[7§, [7§, |83J, gg, [ST], [TJJ, fng, [TTg [12^, TO TO P3j, g3J, f3_7|, 



T481, EH, EH> EH EH EH EH EH' EH EH B> H EH £H> HH EH' 



215i HI HI, |20|, [227[ |228|, H [24^, H, |42|, |24|, |4|, |47|, H |250], |55 



p5B] |257| , |5|, [259], g64[ |272|, gTpJ, [274], |79|, |78|, |87|, |8|, ^89), |91|, |29|, |93|, 
|296| , |9l, |29^, |99l |U|, |03|, |U|, gl5|, |3T|, gl^, |T|, ^22j, ggg, |324[, [~. 
|32G| |327| , p30| , pi], p5] |336[ , pTj [338| , [339j |340[ , |341] [342|, [343] p4^, |352|, [~ ■ 
p57j |359| , g6C| , |36ll , [362| |364] , |365l [366| , |367j [368| , |369l f370| , ^72j |7|, |75|, |97|, 
409] p^ , p9| , p4| . 



3.2 Historical Period 

Of course the division into pre-1981 and post-1981 articles is at a deeper level 
somewhat deceptive. There have been several analogue models investigated 
over the years, with different aims, different levels of sophistication, and 
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ultimately different levels of development. Armed with a good library and 
some hindsight it is possible to find interesting analogues in a number of 
places. 17 

3.2.1 Optics 

Perhaps the first paper to seriously discuss analogue models and effective 
metric techniques was that of Gordon (yes, he of the Klein-Gordon equa- 
tion) [|1 53|| Note that Gordon seemed largely interested in trying to describe 



dielectric media by an "effective metric". That is: Gordon wanted to use a 
gravitational field to mimic a dielectric medium. What is now often referred 
to as the Gordon metric is the expression 

[^effective] fj,v = + i 1 ~ U ~ 2 ] V /jVv> ( 116 ) 

where r]^ u is the flat Minkowski metric, n is the refractive index, and is 
the 4- velocity of the medium. 

After that, there was sporadic interest in effective metric techniques. One 
historically important contribution was one of the problems in the well-known 
book "The classical theory of fields" by Landau and Lifshitz | j224j| . See the 
end of chapter 10, paragraph 90, and the problem immediately thereafter: 
"Equations of electrodynamics in the presence of a gravitational field" . Note 
that in contrast to Gordon, here the interest is in using dielectric media to 
mimic a gravitational field. 

In France the idea was taken up by Pham Mau Quan ||310|| , who showed 



that (under certain conditions) Maxwell's equations can be expressed directly 
in terms of the effective metric specified by the coefficients 



[^effective 



1 

1 

e/ij 



V,V U , (117) 



where g^ v is the ordinary spacetime metric, e and fi are the permeability and 
permittivity, and is the 4-velocity of the medium. The trajectories of the 

17 Indeed historically, though not of direct relevance to general relativity, analogue mod- 
els played a key role in the development of clcctromagnctism — Maxwell's derivation of 
his equations for the electromagnetic field was guided by a rather complicated analogue 
model in terms of spinning vortices of aether. Of course, once you have the equations in 
hand you can treat them in their own right and forget the model that guided you — which 
is exactly what happened in this particular case. 
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electromagnetic rays are interpreted in this case as geodesies of null length 
of this new effective metric. 

Three articles that directly used the dielectric analogy to analyze spe- 
cific physics problems are those of Skrotskii ||354|| , Balazs ||, and Winter- 



berg |J 8[. The general formalism was more fully developed in articles such 



as those by Peblanski [306, 305 , and good summary of this classical period 



can be found in the article by de Felice ||118| . 

In summary and with the benefit of hindsight: An arbitrary gravitational 
field can always be represented as an equivalent optical medium, but subject 
to the somewhat unphysical restriction that 

[magnetic permitivity] oc [electric permeability]. (H8) 

If an optical medium does not satisfy this constraint (with a position inde- 
pendent proportionality constant) then it is not completely equivalent to a 
gravitational field. For a position dependent proportionality constant com- 
plete equivalence can be established in the geometric optics limit, but for 
wave optics the equivalence is not complete. 

3.2.2 Acoustics 

There were several papers in the 1980's using an acoustic analogy to inves- 
tigate the propagation of Shockwaves in astrophysical situations, most no- 
tably those of Moncrief |7(J and Matarrese gkj |262], [26TJ. (Logically they 



are orthogonal to modern developments and belong to the historical period. 
Temporally several of these articles post-date Unruh's 1981 paper by a few 
years, but there seems to have not been any cross-connection.) 

3.2.3 Electro-mechanical analogy 

The so-called "electro-mechanical analogy" has also had a long history within 
the engineering community. It is sometimes extended to obtain an "electro- 
mechanical-acoustic" analogy, or even an "electro-thermal" analogy. Unfor- 
tunately the issues of interest to the engineering community rarely resonate 
within the relativity community, and these engineering analogies (though 
powerful in their own right) have no immediate impact for our purposes. 18 



18 



A recent attempt at connecting the electro-mechanical analogy back to relativity can 



be found in [ 4 3 5 1 
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3.3 Motivation 

The motivation for these investigations (both historical and current) is rather 
mixed. In modern language the reasons to investigate analogue models are: 

• Partly to use condensed matter to gain insight into classical general 
relativity. 

• Partly to use condensed matter to gain insight into curved-space quan- 
tum field theory. 

• Partly to develop an observational window on curved-space quantum 
field theory. 

• Partly to use classical general relativity to gain insight into condensed 
matter physics. 

• Partly (much more tenuous) to gain insight into new and radically 
different ways of dealing with "quantum gravity" . 

3.4 Going further 

There is not much more that we can usefully say here. We have doubtless 
missed some articles of historical importance, but with a good library or a 
fast internet connection the reader will be in as good a position as we are to 
find any additional historical articles. 
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4 A catalogue of models 



In this chapter we will attempt to categorize the very many analogue models 
researchers have investigated. Perhaps the most basic subdivision is into 
classical models and quantum models, but even then many other levels of 
refinement are possible. Consider for instance the following list: 

• Classical models 

— Classical sound. 

— Water waves (gravity waves). 

— Classical refractive index. 

— Normal modes. 

• Quantum models: 

— Bose-Einstein condensates [BECs]. 

— The Heliocentric universe. 

(Helium as an exemplar for just about anything.) 

— Slow light. 

We will now provide a few words on each of these topics. 

4.1 Classical models 
4.1.1 Classical sound 

Sound in a moving fluid has already been extensively discussed in section 
(0), and we will not repeat such discussion here. In contrast, sound in a solid 
exhibits its own distinct and interesting features, notably in the existence of 
a generalization of the normal notion of birefringence — longitudinal modes 
travel at a different speed (typically faster) than do transverse modes. This 
may be viewed as an example of an analogue model which breaks the "light 
cone" into two at the classical level; as such this model is not particularly 
useful if one is trying to simulate special relativistic kinematics with its uni- 
versal speed of light, though it may be used to gain insight into yet another 
way of "breaking" Lorentz invariance. 
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4.1.2 Shallow water waves (gravity waves) 



A wonderful example of the occurrence of an effective metric in nature is 
that provided by gravity waves in a shallow basin filled with liquid 7|. 19 



(See figure |7|). If one neglects the viscosity and considers an irrotational flow, 
v = V0, one can write Bernoulli's equation in the presence of Earth's gravity 

as 



^ + ^(V0) 2 = -^-^-V]|. (119) 

Here p is the density of the fluid, p its pressure, g the gravitational accel- 
eration and VTi a potential associated with some external force necessary to 

establish an horizontal flow in the fluid. We denote that flow by v^. We must 
also impose the boundary conditions that the pressure at the surface, and the 
vertical velocity at the bottom, both vanish: p(z = hs) = 0, v±(z = 0) = 0. 

Once a horizontal background flow is established, one can see that the 
perturbations of the velocity potential satisfy 

dtSt + vi- V||<ty = ~. (120) 
If we now expand this perturbation potential in a Taylor series 

6(f> = y2^5 ( f> n (x,y), (121) 

n=0 



it is not difficult to prove J347J] that surface waves with long wavelengths, 



(long compared with the depth of the basin, A 3> /is), can be described to a 
good approximation by S(po(x, y) and that this field "sees" an effective metric 
of the form 



ds 2 = — 



(c 2 - v)j) dt 2 - 2w\ ■ dx dt + dx ■ dx 



122) 



19 Of course we now mean "gravity wave" in the traditional fluid mechanics sense of 
a water wave whose restoring force is given by ordinary Newtonian gravity. Waves in 
the fabric of spacetime are more properly called "gravitational waves" , though this usage 
seems to be in decline within the general relativity community. Be very careful in any 
situation where there is even a possibility of confusing the two concepts. 
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where c = y 'ghs- The link between small variations of the potential field and 
small variations of the position of the surface is provided by the following 
equation 

5 V± = -h B V|<ty = d t 5h + v| ■ V l{ 6h = ^- Sh. (123) 

The entire previous analysis can be generalized to the case in which the 
bottom of the basin is not flat, and the background flow not purely horizontal 



347 1 . Therefore, one can create different effective metrics for gravity waves 
in a shallow fluid basin by changing (from point to point) the background 
flow velocity and the depth, hs(x, y). 




Figure 7: Gravity waves in a shallow fluid basin with a background horizontal 
flow. 

The main advantage of this model is that the velocity of the surface waves 
can very easily be modified by changing the depth of the basin. This velocity 
can be made very slow, and consequently, the creation of ergoregions should 
be relatively easier than in other models. As described here, this model 
is completely classical, so cannot serve us to study Hawking evaporation an 
other quantum phenomena. However, it can certainly serve to investigate the 
classical phenomena of mode mixing that underlies the quantum processes. 
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4.1.3 Classical refractive index 

The macroscopic Maxwell equations inside a dielectric take the well-known 
form 

V ■ B = 0, V x E + d t B = 0, (124) 

V ■ D = 0, V x H - d t I> = 0, (125) 

with the constitutive relations H = ax 1 ■ B and D = e ■ E. Here, e is the 
3x3 permittivity tensor and n the 3x3 permeability tensor of the medium. 
These equations can be written in a condensed way as 

8 a (Z^ F vP ) = (126) 

where F u p = A{ Vj @\ is the electromagnetic tensor, 

Fa = —F i0 = —Ei, Fij = EijkB k , (127) 

and (assuming the medium is at rest) the non-vanishing components of the 
4'th rank tensor Z are given by 

gOiOj _ _2?0ij0 _ giOjO — _£iO0j _ _}_ e ij. (128) 

Z ijkl = _ £ ijm £ kln ^-4 . ( 12 g) 

supplemented by the conditions that Z is antisymmetric on its first pair of 
indices and antisymmetric on its second pair of indices. Without significant 
loss of generality we can ask that Z also be symmetric under pairwise inter- 
change of the first pair of indices with the second pair — thus Z exhibits 
most of the algebraic symmetries of the Riemann tensor, though this appears 
to merely be accidental, and not fundamental in any way. 

If we compare this to the Lagrangian for electromagnetism in curved 
spacetime 

£ = V^g g m g vP F af3 (130) 

we see that in curved spacetime we can also write the electromagnetic equa- 
tions of motion in the form ( |1 26| ) where now (for some constant K): 
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If we consider a static gravitational field we can always re-write it as a con- 
formal factor multiplying an ultra-static metric 



= ^ {-I ® gij} (132) 

then 

z oiQj = _ z otjo = Z tOjo = _ z mj = _ R g ij. ( 133 ) 

Z^ 1 = K a/— g {g ik g ]l - g il g jk } (134) 

The fact that Z is independent of the conformal factor Q is simply the 
reflection of the well-known fact that the Maxwell equations are conformally 
invariant in (3+1) dimensions. Thus if we wish to have the analogy (be- 
tween a static gravitational field and a dielectric medium at rest) hold at the 
level of the wave equation (physical optics) we must satisfy the two stringent 
constraints 

K gU = l - e«; (135) 

K {g lk g Jl - g u g ]k ) = \ e ijm e kln (136) 

The second of these constraints can be written as 

K V=g s ijm sun {g ik g jl } = ACn- (137) 
In view of the standard formula for 3 x 3 determinants 

Eijm e k in {X %k X? 1 } = 2 detX X~l (138) 

this now implies 



= /i" 1 , (139) 



whence 



Comparing this with 



2^ v 7 ^ g ij = (140) 



2K yf^ J* = (141) 
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we now have: 



,'J 



■i'J 



4 K 2 
AK 2 ,, 



dete 



4 K 2 det/x 



(142) 
(143) 

(144) 



To rearrange this, introduce the matrix square root [/x 1//2 ] lJ , which always 
exists because fj, is real positive definite and symmetric. Then 



^ 2 6 ^' 2 1 V2 1 V 

det(/x e) 



(145) 



Note that if you are given the static gravitational field (in the form Q, gij) 
one can always solve to find an equivalent analogue in terms of permittiv- 
ity/permeability (albeit an analogue that satisfies the mildly unphysical con- 
straint e oc /J.). 20 On the other hand, if you are given permeability and 
permittivity tensors e and n, then it is only for that subclass of media that 
satisfy e oc /j, that one can perfectly mimic all of the electromagnetic effects 
by an equivalent gravitational field. Of course this can be done provided 
one only considers wavelengths that are sufficiently long for the macroscopic 
description of the medium to be valid. In this respect it is interesting to note 
that the behaviour of the refractive medium at high frequencies has been 
used to introduce an effective cutoff for the modes involved in Hawking ra- 
diation |~ H . We shall encounter this model later on when we shall consider 



the trans-Planckian problem for Hawking radiation. 



Eikonal approximation: With a bit more work this discussion can be ex- 
tended to a medium in motion, leading to an extension of the Gordon metric. 
Alternatively, one can agree to ask more limited questions by working at the 
level of geometrical optics (adopting the eikonal approximation), in which 
case there is no longer any restriction on the permeability and permittivity 
tensors. To see this construct the matrix 

20 The existence of this constraint has been independently re-derived several times in 
the literature. In contrast, other segments of the literature seem blithely unaware of this 
important restriction on just when permittivity and permeability arc truly equivalent to 
an effective metric. 
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(146) 



The dispersion relations for the propagation of photons (and therefore the 
sought for geometrical properties) can be obtained from the reduced deter- 
minant of C (notice that the [full] determinant of C is identically zero as 
C fiu k l/ = 0; the reduced determinant is that associated with the three direc- 
tions orthogonal to k v = 0). By choosing the gauge A = one can see that 
this reduced determinant can be obtained from the determinant of the 3x3 
sub- matrix C l K This determinant is 

det(C^) = ^det {-u 2 e ij + e ikm e 3ln ^ n k k h) . (147) 
8 

or, after making some manipulations, 

det(C iJ ) = ^det \-uo 2 e i] + (det/*) -1 ^ n kl k k h - ^k^h)] . (148) 
8 

To simplify this, again introduce the matrix square roots [fj, 1 ^ 2 ]' 13 and 
[/t -1 / 2 ]^-, which always exist because the relevant matrices are real positive 
definite and symmetric. Then define 

k { = {ti 1/2 ] lj kj (149) 

and 

[e] y = det(Ai) [/*- 1/2 6M" 1/2 ]y (150) 

so that 

det(C^) oc det j-cj 2 [e] lj + 5 lj [S mn ~k m k n ] - k l k j }. (151) 

The behaviour of this dispersion relation now depends critically on the 
way that the eigenvalues of e are distributed. 

3 degenerate eigenvalues: If all eigenvalues are degenerate then e = 
e I, implying e oc fi but now with the possibility of a position dependent 
proportionality factor (in the case of physical optics the proportionality factor 
was constrained to be a position-independent constant). In this case we now 
easily evaluate 
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while 

That is 
with 



tr(e) _ tr(e) 

e= — ; — - u and e = detzz — ; — -, (152) 

tr(/x)^ ^tr(/x)' 1 ; 

det(C ii ) oc uo 2 {uo 2 - [r l S mn k m ~k n }} 2 . (153) 



det(C^) oc uo 2 {uo 2 - [g ij k { kj}} 2 , (154) 

u = 1 , }ij = tr(M) W = tr(e) W (155) 
e L J tr(e) det/x tr(/x) dete' v ; 

This last result is compatible with but more general than the result ob- 
tained under the more restrictive conditions of physical optics. In the situa- 
tion where both permittivity and permeability are isotropic, (e iJ = e and 

= fi 5^) this reduces to the perhaps more expected result 

9 lJ = — • (156 
e [i 

2 degenerate eigenvalues: If e has two distinct eigenvalues then the de- 
terminant det(C u ) factorizes into a trivial factor of uo 2 and two quadratics. 
Each quadratic corresponds to a distinct effective metric. This is the physical 
situation encountered in uni-axial crystals, where the ordinary and extraor- 
dinary rays each obey distinct quadratic dispersion relations |39|. From the 
point of view of analogue models this corresponds to a two-metric theory. 

3 distinct eigenvalues: If e has three distinct eigenvalues then the deter- 
minant det(C*- 7 ) is the product of a trivial factor of uj 2 and a non-factorizable 
quartic. This is the physical situation encountered in bi-axial crystals |39|, 



400 , and it seems that no meaningful notion of effective Riemannian metric 



can be assigned to this case. (The use of Finsler geometries in this situation 



is an avenue that may be worth pursuing [. 



Abstract linear electrodynamics: Hehl and co-workers have champi- 
oned the idea of using the linear constitutive relations of electrodynamics as 
the primary quantities, and then treating the spacetime metric (even for flat 
space) as a derived concept. See ||290| , |169| , |223| , |170| . 
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Nonlinear electrodynamics: In general, the permittivity and permeabil- 
ity tensors can be modified by applying strong electromagnetic fields (this 
produces an effectively non-linear electrodynamics). The entire previous dis- 
cussion still applies if one considers the photon as the linear perturbation of 
the electromagnetic field over a background configuration 



(157) 



The background field F}* sets the value of e u (F bg ), and /j,^(F hg ). Equation 
( |126| ) then becomes an equation for fff. This approach has been extensively 
investigated by Novello and co-workers pM |85L pM |8|, |82L |8l ff2| . 



Summary: The propagation of photons in a dielectric medium character- 
ized by 3 x 3 permeability and permittivity tensors constrained by e oc fi is 
equivalent (at the level of geometric optics) to the propagation of photons in 
a curved spacetime manifold characterized by the ultrastatic metric ( |155| ), 
provided one only considers wavelengths that are sufficiently long for the 
macroscopic description of the medium to be valid. If in addition one takes a 
fluid dielectric, by controlling its flow one can generalize the Gordon metric 
and again reproduce metrics of the Painleve-Gullstrand type, and therefore 
geometries with ergoregions. If the proportionality constant relating e oc /x 
is position independent, one can make the stronger statement ( 144 ) which 
holds true at the level of physical optics. 



4.1.4 Normal mode meta-models 

We have already seen how linearizing the Euler-Lagrange equations for a sin- 
gle scalar field naturally leads to the notion of an effective spacetime metric. 
If more than one field is involved the situation becomes more complicated, 
in a manner similar to that of geometrical optics in uni-axial and bi-axial 
crystals. (This should, with hindsight, not be too surprising since electro- 
magnetism, even in the presence of a medium, is definitely a Lagrangian 
system and definitely involves more than one single scalar field.) A normal 
mode analysis based on a general Lagrangian (many fields but still first order 
in derivatives of those fields) leads to a concept of refringence, or more specif- 
ically multi-refringence, a generalization of the birefringence of geometrical 
optics. To see how this comes about, consider a straightforward generaliza- 
tion of the one-field case. 
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We want to consider linearized fluctuations around some background so- 
lution of the equations of motion. As in the single-field case we write (here 
we will follow the notation and conventions of ftTHI) 



<f) A (t, x) = cj> A (t, x) + e cj> A (t, x) + - <f) A (t, x) + 0(e 3 ). 



Now use this to expand the Lagrangian 



lA iA\ 



+ 



+ 



_ lA dC lA 



dC Q J.A , 9C - 



+ 



d 2 C 

ld(d,<j) A ) d(d„<f> B ) 
d 2 C 



U ^2 



d^f d v <t>i 



B 



d(d^ A ) dct> 
+0(e 3 ). 



a, 



158) 



(159) 



Consider the action 



& d+1 xC(d,<t> A A A )- 



(160) 



Doing so allows us to integrate by parts. As in the single-field case we can 
use the Euler-Lagrange equations to discard the linear terms (since we are 
linearizing around a solution of the equations of motion) and so get 



S\ 



id+1 



X 



d 2 c 



d{d^ A ) d(d u <p B 



d^f d v <j> B 
d 2 C 



0(e 



(161) 



Because the fields now carry indices (AB) we cannot cast the action into 
quite as simple a form as was possible in the single-field case. The equation 
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of motion for the linearized fluctuations are now read off as 



- d ^ wS^ ~(ww)^=°- (162) 

This is a linear second-order system of partial differential equations with 
position-dependent coefficients. This system of PDEs is automatically self- 
adjoint (with respect to the trivial "flat" measure d d+1 x). 

To simplify the notation we introduce a number of definitions. First 

f / d 2 C d 2 C \ 

rAB = 2 \d(d^ A ) 0(8^) + djd^) W~^)J ■ (163) 

This quantity is independently symmetric under interchange of /x, v and A, B. 
We will want to interpret this as a generalization of the "densitized metric" , 
f^ v ', but the interpretation is not as straightforward as for the single-field 
case. Next, define 



P' 



d 2 C d 2 C 



AB 



8(8^) d(j) B d{d ll <p B ) d<f) A 

+ U ( d l£ ^ (164 ) 

^ 2 " \d(d u cj> A ) d(d^ B ) d(d^ A ) d(d u 4> B )J ' 1 ; 

This quantity is anti-symmetric in A, B. One might want to interpret this 
as some sort of "spin connexion", or possibly as some generalization of the 
notion of "Dirac matrices" . Finally, define 

d 2 C 1 / d 2 C \ I ( d 2 C \ 

Kab = ~d<f> A d<f> B + 2^ U(«9^ A ) d<p B ) + 2^ \d(d,<P B ) d<p A ) ' (165) 

This quantity is by construction symmetric in (AB). We will want to in- 
terpret this as some sort of "potential" or "mass matrix". Then the crucial 
point for the following discussion is to realize that equation ( | . : |) can be 
written in the compact form 

dp (F u ab d v <j> B ) + i [Y% d^f + d,(Y% 0f )] + K AB <j> B = 0. (f 66) 
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Now it is more transparent that this is a formally self-adjoint second-order 
linear system of PDEs. Similar considerations can be applied to the lineariza- 
tion of any hyperbolic system of second-order PDEs. 

Consider an eikonal approximation for an arbitrary direction in field 
space, that is, take 

(j) A (x) = e A (x) exp[-z^(s)], (167) 

with e A (x) a slowly varying amplitude, and <p(x) a rapidly varying phase. 
In this eikonal approximation (where we neglect gradients in the amplitude, 
and gradients in the coefficients of the PDEs, retaining only the gradients of 
the phase) the linearized system of PDEs (|166|) becomes 



{F u ab d^{x) d v <p{x) + T>1 ab d^{x) + K AB } e B {x) = 0. (168) 

This has a nontrivial solution if and only if e A (x) is a null eigenvector of the 
matrix 

F v ab K K + T» AB k^ + K AB , (169) 

where = d^ip^x). Now, the condition for such a null eigenvector to exist 
is that 

F(x, k) = det {r AB k v + T» AB k^ + K AB ) = 0, (170) 

with the determinant to be taken on the field space indices AB. This is 
the natural generalization to the current situation of the Fresnel equation 
of bi-refringent optics p9| , |226|| . Following the analogy with the situation 



in electrodynamics (either nonlinear electrodynamics, or more prosaically 
propagation in a bi-refringent crystal), the null eigenvector e A (x) would cor- 
respond to a specific "polarization". The Fresnel equation then describes 
how different polarizations can propagate at different velocities (or in more 
geometrical language, can see different metric structures). In particle physics 
language this determinant condition F(x, k) = is the natural generalization 
of the "mass shell" constraint. Indeed it is useful to define the mass shell as 
a subset of the cotangent space by 

F{x) = [k^ I F(x,k) = o}. (171) 

In more mathematical language we are looking at the null space of the de- 
terminant of the "symbol" of the system of PDEs. By investigating F(x, k) 
one can recover part (not all) of the information encoded in the matrices 
f^ u AB , T^ABi an d K AB , or equivalently in the "generalized Fresnel equation" 
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( |170| ). (Note that for the determinant equation to be useful it should be non- 
vacuous; in particular one should carefully eliminate all gauge and spurious 
degrees of freedom before constructing this "generalized Fresnel equation" , 
since otherwise the determinant will be identically zero.) We now want to 
make this analogy with optics more precise, by carefully considering the no- 
tion of characteristics and characteristic surfaces. We will see how to extract 
from the the high-frequency high-momentum regime described by the eikonal 
approximation all the information concerning the causal structure of the the- 
ory. 

One of the key structures that a Lorentzian spacetime metric provides is 
the notion of causal relationships. This suggests that it may be profitable to 
try to work backwards from the causal structure to determine a Lorentzian 
metric. Now the causal structure implicit in the system of second-order PDEs 
given in equation (|166|) is described in terms of the characteristic surfaces, 



and it is for this reason that we now focus on characteristics as a way of encod- 
ing causal structure, and as a surrogate for some notion of Lorentzian metric. 
Note that via the Hadamard theory of surfaces of discontinuity the charac- 
teristics can be identified with the infinite-momentum limit of the eikonal 



approximation [157]. That is, when extracting the characteristic surfaces 
we neglect subdominant terms in the generalized Fresnel equation and focus 
only on the leading term in the symbol (f^ v ab)- in particle physics language 
going to the infinite-momentum limit puts us on the light cone instead of the 
mass shell; and it is the light cone that is more useful in determining causal 
structure. The "normal cone" at some specified point x, consisting of the 
locus of normals to the characteristic surfaces, is defined by 

AT(x) = {fc„ | det {r AB k, fc M ) = 0} . (172) 

As was the case for the Fresnel equation ( |170| ), the determinant is to 
be taken on the field indices AB. (Remember to eliminate spurious and 
gauge degrees of freedom so that this determinant is not identically zero.) 
We emphasise that the algebraic equation defining the normal cone is the 
leading term in the Fresnel equation encountered in discussing the eikonal 
approximation. If there are N fields in total then this "normal cone" will 
generically consist of N nested sheets each with the topology (not necessarily 
the geometry) of a cone. Often several of these cones will coincide, which is 
not particularly troublesome, but unfortunately it is also common for some 
of these cones to be degenerate, which is more problematic. 
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It is convenient to define a function Q(x, k) on the co-tangent bundle 



Q(x,k)=det(f^ AB (x) k, k N 



(173) 



The function Q(x, k) defines a completely-symmetric spacetime tensor (ac- 
tually, a tensor density) with 2N indices 



Q{x,k) = Q W 



filial— Kivfjv | 



*^) ^/xi k^ k V2 



[in) 



(Remember that ab is symmetric in both /ii/ and AB independently.) 
Explicitly, using the expansion of the determinant in terms of completely 
antisymmetric field-space Levi-Civita tensors 



Q 



e A 1 A 2 A 3 -A N B 1 B 2 B 3 -B N 



;i75) 



J J A 2 B 2 J A N B h 



In terms of this Q(x, k) function the normal cone is 



N(x) 



Q{x,k) = 



176) 



In contrast, the "Monge cone" (aka "ray cone", aka "characteristic cone", 
aka "null cone" ) is the envelope of the set of characteristic surfaces through 
the point x. Thus the "Monge cone" is dual to the "normal cone" , its explicit 
construction is given by (Courant and Hilbert |Jjl, volume 2, pp 583): 



Mix) 



t" 



dQ(x, k) 
dk„. 



k^N{ 



x 



(177) 



The structure of the normal and Monge cones encode all the information 
related with the causal propagation of signals associated with the system of 
PDEs. We will now see how to relate this causal structure with the existence 
of effective spacetime metrics, from the experimentally favoured single-metric 
theory compatible with the Einstein equivalence principle to the most com- 
plicated case of pseudo-Finsler geometries ||187[ . 



Suppose that J^ab factorizes 



AB 



h AB r. 



Then 



Q(x,k) = det(h AB ) [r K K 



N 



(178) 
(179) 
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The Monge cones and normal cones are then true geometrical cones 
(with the N sheets lying directly on top of one another). The normal 
modes all see the same spacetime metric, defined up to an unspecified 
conformal factor by g^ v oc This situation is the most interesting 
from the point of view of general relativity Physically it corresponds 
to a single-metric theory, and mathematically it corresponds to a strict 
algebraic condition on the }'^ u ab- 

• The next most useful situation corresponds to the commutativity con- 
dition: 

r ab r p bc = r p ab r bc\ that is [/^,/ a/J ] = o. (iso) 

If this algebraic condition is satisfied, then for all spacetime indices \iv 
and a/3 the ab can be simultaneously diagonalized in field space 
leading to 

Tab = diag{/r, iT, jf, • • • /H (181) 

and then 

JV 

Q(x,k) = H ifTKK] (182) 

A=l 

This case corresponds to an iV-metric theory, where up to an unspec- 
ified conformal factor g^f oc fjf. This is the natural generalization of 
the two metric situation in bi-axial crystals. 

• If f^AB is completely general, satisfying no special algebraic condi- 
tion, then Q(x, k) does not factorize and is in general a polynomial of 
degree 2N in the wave vector k^. This is the natural generalization 
of the situation in bi-axial crystals. (And for any deeper analysis of 
this situation one will almost certainly need to adopt pseudo-Finsler 
techniques |[L87|1 .) 

The message to be extracted from this rather formal discussion is that 
effective metrics are rather general and mathematically robust objects that 
can arise in quite abstract settings — in the abstract setting discussed here 
it is the algebraic properties of the object f^ v ab that eventually leads to 
mono-metricity, multi-metricity, or worse. The current abstract discussion 
also serves to illustrate, yet again, 
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1. that there is a significant difference between the levels of physical nor- 
mal modes (wave equations), and geometrical normal modes (dispersion 
relations), and 

2. that the densitized inverse metric is in many ways more fundamental 
than the metric itself. 



4.2 Quantum models 

4.2.1 Bose— Einstein condensates 

We have seen that one of the main aims of research in analogue models of 
gravity is the possibility of simulating semiclassical gravity phenomena, such 
as the Hawking radiation effect or cosmological particle production. In this 
sense systems characterized by a high degree of quantum coherence, very 
cold temperatures, and low speeds of sound offer the best test field. Hence 
it is not surprising that in recent years Bose-Einstein condensates (BECs) 
have become the subject of extensive study as possible analogue models of 
general relativity |3§, ggg, g§ gg, gg, |TT6;, [TT§. 

Let us start by very briefly reviewing the derivation of the acoustic met- 
ric for a BEC system, and show that the equations for the phonons of the 
condensate closely mimic the dynamics of a scalar field in a curved space- 
time. In the dilute gas approximation, one can describe a Bose gas through 
a quantum field \l/ satisfying 




■c(a) J (183) 



Here k parameterizes the strength of the interactions between the different 
bosons in the gas. It can be re-expressed in terms of the scattering length as 

. , 4nah 2 

K (a) = . 184 

m 

As usual, the quantum field can be separated into a macroscopic (classical) 
condensate and a fluctuation: \& = ip + tp, with = ip. Then, by adopting 
the self-consistent mean field approximation (see for example ||155|| ) 

(p*(p(p ~ 2{(p i (p) (p + ((pip) <p\ (185) 



(37 



one can arrive at the set of coupled equations: 

ih— ^(t,x) = ( -— V 2 + Kxt(x) + Kn c \ V(t,x) 

+K{2nip(t,x) + rhi>*{t,x)}; (186) 

d ( h 2 \ 

ih—<p(t,x) = f -— V 2 + Kxt(x) + k 2n T J £(t,x) 

+/cmr^(t,x). (187) 

Here 

n c = |^(t,x)| 2 ; m c = ^ 2 (*,x); (188) 

n = {$(p); m={lf(p); (189) 

ny = n c + n; my = m c + rh. (190) 

The equation for the classical wave function of the condensate is closed only 
when the back-reaction effect due to the fluctuations are neglected. (This 
back-reaction is hiding in the parameters rh and n.) This is the approxima- 
tion contemplated by the Gross-Pitaevskii equation. In general one will have 
to solve both equations simultaneously. Adopting the Madelung representa- 
tion for the wave function of the condensate 

V>(t,x) = y/njt^x) exp[-i0(t,x)/fc], (191) 

and defining an irrotational "velocity field" by v = V0/m, the Gross- 
Pitaevskii equation can be rewritten as a continuity equation plus an Euler 
equation: 

^n c + V • (n c v) = 0, (192) 



These equations are completely equivalent to those of an irrotational and 
inviscid fluid apart from the existence of the so-called quantum potential 

^quantum = -?? V 2 y/nj {2my/rQ , (194) 



OS 



which has the dimensions of an energy. Note that 



^i^quantum — f^c 



h 2 V- 



2m 



v.- 



4m 



n c ViV 7 - lnn c 



(195) 



which justifies the introduction of the so-called quantum stress tensor 

h 2 



a 



quantum 



Am 



n c ViVj lnn c . 



(196) 



This tensor has the dimensions of pressure, and may be viewed as an intrin- 
sically quantum anisotropic pressure contributing to the Euler equation. If 
we write the mass density of the Madelung fluid as p = m n c , and use the 
fact that the flow is irrotational then the Euler equation takes the form 



P 



-v+(v.V)v 



p V 



Knrt(t,x) 



m 



+ V 



2m 2 



V . a quantum = Q 



Note that the term V cy± /m has the dimensions of specific enthalpy, while 
Kp 2 / (2m) represents a bulk pressure. When the gradients in the density of 
the condensate are small one can neglect the quantum stress term leading to 
the standard hydrodynamic approximation. Because the flow is irrotational, 
the Euler equation is often more conveniently written in Hamilton-Jacobi 
form: 
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Apart from the wave function of the condensate itself, we also have to account 
for the [typically small] quantum perturbations of the system ( |187|) . These 



quantum perturbations can be described in several different ways, here we 
are interested in the "quantum acoustic representation" 
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where rii,9i are real quantum fields. By using this representation equation 



187|) can be rewritten as 
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Here D 2 represents a second-order differential operator obtained from lin- 
earizing the quantum potential. Explicitly: 
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The equations we have just written can be obtained easily by linearizing 



the Gross-Pitaevskii equation around a classical solution: 



n, 



n c + ni, 



<j) —> (j) + (f>i. It is important to realize that in those equations the back- 
reaction of the quantum fluctuations on the background solution has been 
assumed negligible. We also see in those equations, ( |200| ) and ( |201| ), that 
time variations of V cxt and time variations of the scattering length a appear 
to act in very different ways. Whereas the external potential only influences 
the background equation ( |198|) [and hence the acoustic metric in the analogue 
description], the scattering length directly influences both the perturbation 
and background equations. From the previous equations for the linearized 
perturbations it is possible to derive a wave equation for 9± (or alternatively, 
for ni). All we need is to substitute in equation ( |200| ) the rii obtained from 
equation (|20 1|) . This leads to a PDE that is second-order in time derivatives 
but infinite order in space derivatives — to simplify things we can construct 
the symmetric 4x4 matrix 
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(Greek indices run from 0-3, while Roman indices run from 1-3.) 
introducing (3 + 1) -dimensional space-time coordinates 



(t; x* 



the wave equation for 9\ is easily rewritten as 

d^r dX) = o. 

Where the are differential operators acting on space only: 
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Then, 
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Now, if we make an spectral decomposition of the field Q\ we can see that 
for wavelengths larger than h/mc s (this corresponds to the "healing length", 
as we will explain below), the terms coming from the linearization of the 
quantum potential (the D2) can be neglected in the previous expressions, in 
which case the f^ u can be approximated by numbers, instead of differential 
operators. (This is the heart of the acoustic approximation.) Then, by 
identifying 
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the equation for the field Q\ becomes that of a (massless minimally coupled) 
quantum scalar field over a curved background 
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with an effective metric of the form 

-{c s (a,n c ) 2 -v 2 } 
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Here the magnitude c s (n c ,a) represents the speed of the phonons in the 
medium: 
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With this effective metric now in hand, the analogy is fully established, and 
one is now in a position to start asking more specific physics questions. 



4.2.2 BEC models in the eikonal approximation 

It is interesting to consider the case in which the above "hydrodynamical" 
approximation for BECs does not hold. In order to explore a regime where 
the contribution of the quantum potential cannot be neglected we can use 
the so called eikonal approximation, a high-momentum approximation where 
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the phase fluctuation d\ is itself treated as a slowly-varying amplitude times 
a rapidly varying phase. This phase will be taken to be the same for both ri\ 
and Q\ fluctuations. In fact, if one discards the unphysical possibility that the 
respective phases differ by a time varying quantity, any time-constant differ- 
ence can be safely reabsorbed in the definition of the (complex) amplitudes. 
Specifically, we shall write 
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As a consequence of our starting assumptions, gradients of the amplitude, 
and gradients of the background fields, are systematically ignored relative to 
gradients of <fi. [Warning: what we are doing here is not quite a "standard" 
eikonal approximation, in the sense that it is not applied directly on the 
fluctuations of the field if)(t, x) but separately on their amplitudes and phases 
pi and 0i.] We adopt the notation 
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Then the operator D 2 can be approximated as 
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A similar result holds for D 2 acting on 9\. That is, under the eikonal ap- 
proximation we effectively replace the operator D 2 by the function 
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For the matrix f^ u this effectively results in the replacement 
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(As desired, this has the net effect of making f^ u a matrix of numbers, 
not operators.) The physical wave equation ( | .] ) now becomes a nonlinear 
dispersion relation 



0. 



(225) 



After substituting the approximate D-i into this dispersion relation and re- 
arranging, we see (remember: k 2 = | \k\ \ 2 = 5 lJ ki kj) 
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Introducing the speed of sound c s this takes the form: 



w = «* ki±ilc 2 k 2 +{^k 2 



(228) 



At this stage some observations are in order: 

—(1) It is interesting to recognize that the dispersion relation ( |228| ) is 
exactly in agreement with that found in 1947 by Bogoliubov (reprinted 
in ||312| |; see also ||245|| ) for the collective excitations of a homogeneous Bose 
gas in the limit T — >• (almost complete condensation). In his derivation Bo- 
goliubov applied a diagonalization procedure for the Hamiltonian describing 
the system of bosons. 

—(2) It is easy to see that ( |228| ) actually interpolates between two dif- 
ferent regimes depending on the value of the wavelength A = 27r/||fc|| with 
respect to the "acoustic Compton wavelength" A c = h/(mc s ). (Remember 
that c s is the speed of sound; this is not a standard particle physics Compton 
wavelength.) In particular, if we assume vq = (no background velocity), 
then for large wavelengths A ^> A c one gets a standard phonon dispersion 



73 



relation uo pa c||A;||. For wavelengths A A c the quasi-particle energy tends 
to the kinetic energy of an individual gas particle and in fact u ~ h 2 k 2 / (2m). 

We would like also to highlight that in relative terms, the approximation 
by which one neglects the quartic terms in the dispersion relation gets worse 
as one moves closer to a horizon where = — c s . The non-dimensional 
parameter that provides this information is defined by 



1 + il - 1 1 A 2 

(229) 
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As we will discuss in subsection |5.1.3| , this is the reason why sonic horizons in 



a BEC can exhibit different features from those in standard general relativity. 

—(3) The dispersion relation ( | ; ] ) exhibits a contribution due to the 
background flow v l hi, plus a quartic dispersion at high momenta. The 
group velocity is 

<-fe-4* feM^ (230) 
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Dispersion relations of this type (but in most cases with the sign of the 
quartic term reversed) have been used by Corley and Jacobson in analyzing 
the issue of trans-Planckian modes in the Hawking radiation from general 
relativistic black holes |]188] , |189| , [501. The existence of modified dispersion 
relations (MDR), that is, dispersion relations that break Lorentz invariance, 
can be taken as a manifestation of new physics showing up at high energies/ 
short wavelengths. In their analysis, the group velocity reverses its sign for 
large momenta. (Unruh's analysis of this problem used a slightly different toy 
model in which the dispersion relation saturated at high momentum ||379|| .) 
In our case, however, the group velocity grows without bound allowing high- 
momentum modes to escape from behind the horizon. Thus the acoustic 
horizon is not absolute in these models, but is instead frequency dependent, 
a phenomenon that is common once non-trivial dispersions are included. 

Indeed, with hindsight the fact that the group velocity goes to infinity for 
large k was pre-ordained: After all, we started from the generalized nonlinear 
Schrodinger equation, and we know what its characteristic curves are. Like 
the diffusion equation the characteristic curves of the Schrodinger equation 
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(linear or nonlinear) move at infinite speed. If we then approximate this gen- 
eralized nonlinear Schrodinger equation in any manner, for instance by lin- 
earization, we cannot change the characteristic curves: For any well behaved 
approximation technique, at high frequency and momentum we should re- 
cover the characteristic curves of the system we started with. However, what 
we certainly do see in this analysis is a suitably large region of momentum 
space for which the concept of the effective metric both makes sense, and 
leads to finite propagation speed for medium-frequency oscillations. 

This type of superluminal dispersion relation has also been analyzed by 
Corley and Jacobson p2| . They found that this escape of modes from behind 
the horizon often leads to self-amplified instabilities in systems possessing 
both an inner horizon as well as an outer horizon, possibly causing them to 
disappear in an explosion of phonons. This is also in partial agreement with 
the stability analysis performed by Garay et al | |138| , |139|| using the whole 



Bogoliubov equations. Let us however leave further discussion regarding 



these developments to the subsection [5.1.3| on horizon stability. 



4.2.3 The Heliocentric universe 

Helium is one of the most fascinating elements provided by nature. Its struc- 
tural richness confers on helium a paradigmatic character regarding the emer- 
gence of many and varied macroscopic properties from the microscopic world 
(see |f419|| and references therein). Here, we are interested in the emergence 
of effective geometries in helium, and their potential use in testing aspects 
of semiclassical gravity. 

Helium four, a bosonic system, becomes superfluid at low temperatures 
(2.17 K at vapour pressure). This superfluid behaviour is associated with 
the condensation in the vacuum state of a macroscopically large number of 
atoms. A superfluid is automatically an irrotational and inviscid fluid, so 
in particular one can apply to it the ideas worked out in Section ||. The 
propagation of classical acoustic waves (scalar waves) over a background 
fluid flow can be described in terms of an effective Lorentzian geometry: the 
acoustic geometry. However, in this system one can naturally go considerably 
further, into the quantum domain. For long wavelengths, the quasiparticles in 
this system are quantum phonons. One can separate the classical behaviour 
of a background flow (the effective geometry) from the behaviour of the 
quantum phonons over this background. In this way one can reproduce, in 
laboratory settings, different aspects of quantum field theory over curved 
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backgrounds. The speed of sound in the superfluid phase is typically of 
the order of cm/ sec. Therefore, at least in principle, it should not be too 
difficult to establish configurations with supersonic flows and their associated 
ergoregions. 

Helium three, the fermionic isotope of helium, in contrast becomes su- 
perfluid at very much lower temperatures (below 2.5 milli-K). The reason 
behind this rather different behaviour is the pairing of fermions to form ef- 
fective bosons (Cooper pairing), which are then able to condense. In the 
so-called 3 He-A phase, the structure of the fermionic vacuum is such that it 
possesses two Fermi points, instead of the more typical Fermi surface. In an 
equilibrium configuration one can choose the two Fermi points to be located 
at {Px = 0,p y = 0,p z = ±Pf} (in this way, the z axis signals the direc- 
tion of the angular momentum of the pairs). Close to either Fermi point 
the spectrum of quasiparticles becomes equivalent to that of Weyl fermions. 
From the point of view of the laboratory, the system is not isotropic, it is 
axisymmetric. There is a speed for the propagation of quasiparticles along 
the z axis, cy ~ cm/ sec, and a different speed, c± ~ 10~ 5 cy, for propagation 
perpendicular to the symmetry axis. However, from an internal observer's 
point of view this anisotropy is not "real", but can be made to disappear 
by an appropriate rescaling the coordinates. Therefore, in the equilibrium 
case, we are reproducing the behaviour of Weyl fermions over Minkowski 
spacetime. Additionally, the vacuum can suffer collective excitations. These 
collective excitations will be experienced by the Weyl quasiparticles as the 
introduction of an effective electromagnetic field and a curved Lorentzian 
geometry. The control of the form of this geometry provides the sought for 
gravitational analogy. 

Apart from the standard way to provide a curved geometry based on pro- 
ducing non-trivial flows, there is also the possibility of creating topologically 
non-trivial configurations with a built-in non-trivial geometry. For example, 
it is possible to create a domain-wall configuration |[203| , [202|| (the wall con- 
tains the z axis) such that the transverse velocity c± acquires a profile in 
the perpendicular direction (say along the x axis) with c±_ passing through 
zero at the wall (see figure H). This particular arrangement could be used to 
reproduce a black hole-white hole configuration only if the soliton is set up 
to move with a certain velocity along the x axis. This configuration has the 
advantage than it is dynamically stable, for topological reasons, even when 
some supersonic regions are created. 

A third way in which superfluid Helium can be used to create analogues 



76 



Figure 8: Domain wall configuration in 3 He. 



of gravitational configurations is the study of surface waves (or ripplons) on 
the interface between two different phases of 3 He E71 |418||. In particular, 



if we have a thin layer of 3 He-A in contact with another thin layer of 3 He- 
B, the oscillations of the contact surface "see" an effective metric of the 
form 1(4171, 131 



ds 2 
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(All of this provided that we are looking at wavelengths larger than the layer 
thickness, kh^ <C 1 and fc/ig <C 1.) 
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Figure 9: Ripplons in the interface between two sliding superffuids. 



The advantage of using surface waves instead of bulk waves in superfluids 
is that one could create horizons without reaching supersonic speeds in the 
bulk fluid. This could alleviate the appearance of dynamical instabilities in 
the system, that in this case are controlled by the strength of the interaction 
of the ripplons with bulk degrees of freedom ||417 , [418 . 



4.2.4 Slow light 

The geometrical interpretation of the motion of light in dielectric media leads 
naturally to conjecture that the use of flowing dielectrics might be useful for 
simulating general relativity metrics with ergoregions and black holes. Un- 
fortunately, these types of geometry require flow speeds comparable to the 
group velocity of the light. Since typical refractive indexes in non-dispersive 
media are quite close to unity, it is then clear that it is practically impos- 
sible to use them to simulate such general relativistic phenomena. How- 
ever recent technological advances have radically changed this state of af- 
fairs. In particular the achievement of controlled slowdown of light, down 
to velocities of a few meters per second (or even down to complete rest) 
16 1| , P07| , |52| , P 14] , |311| , [376| , |350|| , has opened a whole new set of possibilities 



regarding the simulation of curved-space metrics via flowing dielectrics. 

But how can light be slowed down to these "snail-like" velocities? The 
key effect used to achieve this takes the name of Electromagnetically Induced 
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transparency (EIT). A laser beam is coupled to the excited levels of some 
atom and used to strongly modify its optical properties. In particular one 
generally choses an atom with two long-lived metastable (or stable) states, 
plus a higher energy state that has some decay channels into these two lower 
states. The coupling of the excited states induced by the laser light can 
affect the transition from a lower energy state to the higher one, and hence 
the capability of the atom to absorb light with the required transition energy. 
The system can then be driven into a state where the transitions between 
each of the lower energy states and the higher energy state exactly cancel out, 
due to quantum interference, at some specific resonant frequency. In this way 
the higher-energy level has null averaged occupation number. This state is 
hence called a "dark state" . EIT is characterized by a transparency window, 
centred around the resonance frequency, where the medium is both almost 
transparent and extremely dispersive (strong dependence on frequency of the 
refractive index). This in turn implies that the group velocity of any light 
probe would be characterized by very low real group velocities (with almost 
vanishing imaginary part) in proximity to the resonant frequency. 

Let us review the most common setup envisaged for this kind of analogue 
model. A more detailed analysis can be found in |[234|| . One can start by 



considering a medium in which an EIT window is opened via some control 
laser beam which is oriented perpendicular to the direction of the flow. One 
then illuminates this medium, now along the flow direction, with some probe 
light (which is hence perpendicular to the control beam). This probe beam 
is usually chosen to be weak with respect to the control beam, so that it 
does not modify the optical properties of the medium. In the case in which 
the optical properties of the medium do not vary significantly over several 
wavelengths of the probe light, one can neglect the polarization and can 
hence describe the propagation of the latter with a simple scalar dispersion 



relation 237, 126 



k 2 - ^ [1 + XH] > (234) 
er 

where x is the susceptibility of the medium, related to the refractive index n 
via the simple relation n = y/1 + \- 

It is easy to see that in this case the group and phase velocities differ 

duo c uj c 

S dk f TT ^ + u L dx' k v/T+X 
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So even for small refractive indexes one can get very low group velocities, 
due to the large dispersion in the transparency window, and in spite of the 
fact that the phase velocity remains very near to c. (The phase velocity is 
exactly c at the resonance frequency Uq). In an ideal EIT regime the probe 
light experiences a vanishing susceptibility \ near the the critical frequency 
Uq, this allows us to express the susceptibility near the critical frequency via 
the expansion 

= _^ (a, _ Uo ) + o [(cu - u f] , (236) 

where a is sometimes called the "group refractive index". The parameter a 
depends on the dipole moments for the transition from the metastable states 
to the high energy one, and most importantly depends on the ratio between 
the probe-light energy per photon, Huj , and the control-light energy per 
atom ||234|| . This might appear paradoxical because it seems to suggest that 



for a dimmer control light the probe light would be further slowed down. 
However this is just an artificial feature due to the extension of the EIT 
regime beyond its range of applicability. In particular in order to be effective 
the EIT requires the control beam energy to dominate any all processes and 
hence it cannot be dimmed at will. 
At resonance we have 

duj c c uj 
v g = — -> — — « -; v ph = - -> c. (237) 
ok 1 + a a k 

We can now generalize the above discussion to the case in which our highly 
dispersive medium flows with a characteristic velocity profile u(x, t). In 
order to find the dispersion relation of the probe light in this case we just 
need to transform the dispersion relation ( |234j) from the comoving frame 
of the medium to the laboratory frame. Let us consider for simplicity a 
monochromatic probe light (more realistically a pulse with a very narrow 
range of frequencies u near uj ). The motion of the dielectric medium creates 
a local Doppler shift of the frequency 

u -> 7 (u - u ■ k) , (238) 

where 7 is the usual relativistic factor. Given that k 2 — uj 2 /c 2 is a Lorentz 
invariant, it is then easy to see that this Doppler detuning affects the disper- 
sion relation ( |234| ) only via the susceptibility dependent term. Given further 
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that in any realistic case one would deal with non-relativistic fluid velocities 
u C ewe can then perform an expansion of the dispersion relation up to 



second order in u/c. Expressing the susceptibility via ( |236|) we can then 
rewrite the dispersion relation in the form ||2371 



knk u — 0, 
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and (most of the relevant articles adopt the signature 
do for this particular section) 
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(240) 
as we also 
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The inverse of this tensor will be the covariant effective metric experienced by 
the probe light, whose rays would then be null geodesies of the line element 
ds 2 = g^dx^dx" . In this sense the probe light will propagate as in a curved 
background. Explicitly one finds the covariant metric to be 
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Several comments are in order concerning the metric ( p42j ). First of all 
it is clear that although more complicated than an acoustic metric it will be 



still possible to cast it into the Arnowitt-Deser-Misner-like form ||393 
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where the effective speed u ef j is proportional to the fluid flow speed u and 
the three space effective metric g e g is (rather differently from the acoustic 
case) non-trivial. 

In any case, the existence of this ADM form already tells us that an 
ergoregion will always appear once the norm of the effective velocity exceeds 
the effective speed of light (which for slow light is approximately c/a where a 
can be extremely large due to the huge dispersion in the transparency window 
around the resonance frequency uq). However a trapped surface (and hence 
an optical black hole) will form only if the inward normal component of the 
effective flow velocity exceeds the group velocity of light. In the slow light 
setup so far considered such velocity turns out to be u = c/(2y/a). 

The realization that ergoregions and event horizons can be simulated via 
slow light may lead one to the (erroneous) conclusion that this is an optimal 
system for simulating particle creation by gravitational fields. However, as 
pointed out by Unruh in [|286| , |38 1|| , such conclusion would turn out to be over- 
enthusiastic. In order to obtain particle creation an inescapable requirement 
is to have so-called "mode mixing", that is, mixing between the positive 
and negative frequency modes of the incoming and outgoing states. This 
is tantamount to saying that there must be regions where the frequency of 
the quanta as seem by a stationary observer at infinity (laboratory frame) 
becomes negative beyond the ergosphere at g 00 = 0. 

In a flowing medium this can in principle occur thanks to the tilting of 
the dispersion relation due to the Doppler effect caused by the velocity of 
the flow equation ( | : ] ) . but this also tell us that the condition u — u ■ k < 
can be satisfied only if the velocity of the medium exceeds |a>o/A;| which is 
the phase velocity of the probe light, not its group velocity. Since the phase 
velocity in the slow light setup we are considering is very close to c, the 
physical speed of light in vacuum, not very much hope is left for realizing 
analogue particle creation in this particular laboratory setting. 

However it was also noticed by Unruh and Schutzhold in ||381|| that a 



different setup for slow light might deal with this and other issues (see \p8'. 
for a detailed summary). In the setup suggested by these authors there are 
two strong background counter-propagating control beams illuminating the 
atoms. The field describing the beat fluctuations of this electromagnetic 
background can be shown to satisfy, once the dielectric medium is in motion, 
the same wave equation as that on a curved background. In this particular 
situation the phase velocity and the group velocity are approximately the 
same, and both can be made small, so that the previously discussed ob- 
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struction to mode mixing is removed. So in this new setup it is concretely 
possible to simulate classical particle creation such as, e.g., super-radiance 
in the presence of ergoregions. 

Nonetheless the same authors showed that this does not open the possi- 
bility for a simulation of quantum particle production (e.g., Hawking radi- 
ation). This is because that effect also requires the commutation relations 
of the field to generate the appropriate zero-point energy fluctuations (the 
vacuum structure) according to the Heisenberg uncertainty principle. This 
is not the case for the effective field describing the beat fluctuations of the 
system we have just described, which is equivalent to saying that it does 
not have a proper vacuum state (i.e. analogue to any physical field). Hence 
one has to conclude that any simulation of quantum particle production is 
precluded. 



4.3 Going further 



We feel that the catalogue we have just presented is reasonably complete and 
covers the key items. For additional background on any of these topics, we 
would suggest sources such as the books "Artificial Black Holes" ||286|| and 
"The universe in a Helium droplet" ||419 [. For more specific detail, check this 



review's bibliography, and use Spires to check for recent developments. 
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5 Lessons from analogue models 



Of course the entire motivation for looking at analogue models is to be able 
to learn more physics. One could start studying analogue models with the 
idea of seeing whether it is possible (either theoretically or in practice) to 
reproduce in the laboratory various gravitational phenomena whose real ex- 
istence in nature cannot be currently checked. These are phenomena that 
surpass our present (and foreseeable) observational capabilities, but yet, we 
believe in their existence because it follows from seemingly strong theoreti- 
cal arguments within the standard frameworks of general relativity and field 
theory in curved space. However, the interest of this approach is not merely 
to reproduce these gravitational phenomena in some formal analogue model, 
but to see which departures from the ideal case show up in real analogue 
models, and to analyse whether similar deviation are likely to appear in real 
gravitational systems. 

When one thinks about emergent gravitational features in condensed mat- 
ter systems, one immediately realizes that these features only appear in the 
low-energy regime of the analogue systems. When the systems are probed at 
high energies (short length scales) the effective geometrical description of the 
analogue models break down, as one starts to be aware that the systems are 
actually composed of discrete pieces (atoms and molecules). This scenario is 
quite similar to what one expects to happen with our geometrical description 
of the Universe, when explored with microscopic detail at the Planck scale. 

That is, the study of analogue models of general relativity is giving us 
insights as to how the standard theoretical picture of different gravitational 
phenomena could change when taking into account additional physical knowl- 
edge coming from the existence of an underlying microphysical structure. 
Quite robustly, these studies are telling us already that the first deviations 
from the standard general relativity picture can be encoded in the form of 
high-energy violations of Lorentz invariance in particle dispersion relations. 
Beyond these first deviations, the analogue models of general relativity pro- 
vide well understood examples (the underlying physics is well known) in 
which a description in terms of fields in curved spacetimes shows up as a 
low-energy-regime emergent phenomena. 

The analogue models are being used to throw light upon these general 
questions through a number of specific routes. Let us now turn to discussing 
several specific physics issues that are being analyzed from this perspective. 
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5.1 Hawking radiation 
5.1.1 Basics 

As is well known, in 1974 Stephen Hawking announced that quantum me- 
chanically even a spherical distribution of matter collapsing to form a black 
hole should emit radiation; with a spectrum approximately that of a black 
body ||164j , |165|| . A black hole will tend to evaporate by emitting particles 



from its horizon toward infinity. Hawking radiation is a quantum-field- in- 
curved-space effect. The existence of radiation emission is a kinematic effect 
that does not rely on Einstein equations. Therefore, one can aim to reproduce 
it in a condensed matter system. Within standard field theory, a minimal 
requirement for having Hawking radiation is the existence in the background 



configuration of an apparent horizon [396]. So, in principle, to be able to 



reproduce Hawking radiation in a laboratory one would need at least two 
requirements: 

• To choose an adequate analogue system; it has to be a quantum ana- 
logue model (see Section [|) such that its description could be separated 
into a classical effective background spacetime plus some standard rel- 
ativistic quantum fields leaving on it (it can happen that the quantum 
fields do not satisfy the appropriate commutation or anti-commutation 
relations ||381|| ). 



• To configure the analogue geometry such that it includes an apparent 
horizon. That is, within an appropriate quantum analogue model, the 
formation of an apparent horizon for the propagation of the quantum 
fields should excite the fields as to result in the emission of a thermal 
distribution of field particles. 

This is a straight and quite naive translation of the standard Hawking 
effect derivation to the condensed matter realm. However, in reality, this 
translation process has to take into account additional issues that, as we are 
trying to convey, instead of problems, are where the interesting physics lies. 

- (1) The effective description of the quantum analogue systems as fields 
over a background geometry breaks down when probed at sufficiently short 
length scales. This could badly influence the main features of Hawking ra- 
diation. In fact, immediately after the inception of the idea that black holes 
radiate, it was realized that there was a potential problem with the calcula- 



tion I 7| . It strongly relies on the validity of quantum field theory on curved 
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backgrounds up to arbitrary high energies. Following a wave packet with a 
certain frequency at future infinity backwards in time, we can see that it had 
to contain arbitrarily large frequency components with respect to a local free 
fall observer (well beyond the Planck scale) when it was close to the horizon. 
In principle any unknown physics at the Planck scale could strongly influ- 
ence the Hawking process so that one should view it with suspicion. This 
is the so-called trans-Planckian problem of Hawking radiation. To create an 
analogue model exhibiting Hawking radiation will be, therefore, equivalent 
to give a solution to the trans-Planckian problem. 

- (2) In order to clearly observe Hawking radiation, one should first be 
sure that there is no other source of instabilities in the system that could mask 
the effect. In analogue models such as liquid Helium or BECs the interaction 
of a radial flow (with speed of the order of the critical Landau speed, that 
in these cases coincides with the sound speed |216| ) with the surface of the 
container (an electromagnetic potential in the BECs case) might cause the 
production of rotons and quantized vortices, respectively. Thus, in order to 
produce an analogue model of Hawking radiation, one has to be somewhat 
ingenious. For example in the liquid Helium case, instead of taking acoustic 
waves in a supersonic flow stream as the analogue model, it is preferable 
to use as analogue model ripplons in the interface between two different 
phases, A and B phases, of Helium-3 ||417|| . Another option is to start from 
a moving domain wall configuration. Then, the topological stability of the 
configuration prevents its destruction when creating an horizon [ [20 2j , [203]| . In 
the case of BECs a way to suppress the formation of quantized vortices is to 
take effectively one-dimensional configurations. If the transverse dimension 
of the flow is smaller than the healing length then there is no space for 



the existence of a vortex |]19| . In either liquid Helium or BECs, there is 
also the possibility of creating an apparent horizon by rapidly approaching a 
critical velocity profile (see figure |II]), but without actually crossing into the 
supersonic regime ||13|| , softening in this way the appearance of dynamical 
instabilities. 

- (3) Real analogue models cannot, strictly speaking, reproduce eternal 
black-hole configurations. An analogue model of a black hole has always to 
be created at some finite laboratory time. Therefore, one is forced to care- 
fully analyze the creation process, as it can greatly influence Hawking effect. 
Depending on the procedure of creation, one could end up in quite different 
quantum states for the field and only some of them might exhibit Hawking 
radiation. This becomes more important when considering that the analogue 
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Figure 10: Velocity profile for a left going flow; the profile is dynamically 
modified with time so that it reaches the profile with a sonic point at the 
asymptotic future. 



models incorporate modified dispersion relations. An inappropriate prepa- 
ration, together with modified dispersion relation effects, could completely 
eliminate Hawking radiation ||382| . 



- (4) Another important issue is the need to characterize "how quantum" 
a specific analogue model is. Even though, strictly speaking, one could say 
that any system undergoes quantum fluctuations, the point is how important 
they are in its description. In trying to build an analogue model of Hawking's 
quantum effect, the relative value of Hawking temperature with respect to 
the environment is going to tell us whether the system can be really thought 
as a quantum analogue model or as effectively classical. For example, in our 
standard cosmological scenario, for a black hole to radiate at temperatures 
higher than that of the Cosmic Microwave Background, ~ 3 K, the black hole 
should have a diameter of the order of micrometers or less. We would have 
to say that such black holes are no longer classical, but semiclassical. The 
black holes for which we have some observational evidence are of much higher 
mass and size, so their behaviour can be thought of as completely classical. 
Estimates of the Hawking temperature reachable in BECs yield T ~ lOOnK 
p9| . This has the same order of magnitude of the temperature of the BECs 
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themselves. This is telling us that regarding the Hawking process, BECs can 
be considered to be highly-quantum analogue models. 

Because of its importance, let us now review what we know about the 
effects of high-energy dispersion relations on the Hawking process. 



5.1.2 Trans-Planckian problem 

We saw in the introduction to this section that the trans-Planckian problem 
of Hawking radiation was one of the strongest motivations for the modern 
research into analogue models of gravity. In fact it was soon realized that such 
models could provide a physical framework within which a viable solution of 
the problem could be found. Let us explain why and how. 

As we have said, the requirement of a reservoir of ultra-high frequency 
modes nearby the horizon seems to indicate a possible (and worrisome) sen- 
sitivity of the black hole radiation to the microscopic structure of spacetime. 
Actually by assuming exact Lorentz invariance one could in principle always 
locally transform the problematic ultra high frequency modes to low energy 
ones via some appropriate Lorentz transformation ||188|| . However in doing 



so it would have to rely on the physics of reference frames moving ultra fast 
with respect to us, as the reference frame needed would move arbitrarily close 
to the speed of light. Hence we would have to apply Lorentz invariance in a 
regime of arbitrary large boots yet untested and in principle never completely 
testable given the non-compactness of the boost subgroup. The assumption 
of an exact boost symmetry is linked to a scale-free nature of spacetime given 
that unbounded boosts expose ultra-short distances. Hence the assumption 
of exact Lorentz invariance needs, in the end, to rely on some idea on the 
nature of spacetime at ultra-short distances. 

It was this type of reasoning that led in the nineties to a careful recon- 
sideration of the crucial ingredients required for the derivation of Hawking 



radiation | 38], |189|, | jj. In particular investigators explored the possibility 
that spacetime microphysics could provide a short distance, Lorentz-breaking 
cutoff, but at the same time leave Hawking's results unaffected at energy 
scales well below that set by the cutoff. 

Of course ideas about a possible cutoff imposed by the discreteness of 
spacetime at the Planck scale had already been discussed in the literature 
well before Unruh's seminal paper ||378|| . However such ideas were running 



into serious difficulties given that a naive short distance cutoff posed on the 
available modes of a free field theory results in a complete removal of the 
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evaporation process (see e.g., [J | and references therein). Indeed there are 
alternative ways through which the effect of the short scales physics could be 
taken into account, and analogue models provide a physical framework where 
these ideas could be put to the test. In fact analogue models provide explicit 
examples of emergent spacetime symmetries, they can be used to simulate 
black hole backgrounds, they may be endowed with quantizable perturbations 
and, in most of the cases, they have a well known microscopic structure. 
Given that Hawking radiation can be, at least in principle, simulated in such 
systems one might ask how and if the trans-Planckian problem is resolved in 
these cases. 

Modified dispersion relations: The general feature that most of the 
work on this subject has focussed on is that in analogue models the quasi- 
particles propagating on the effective geometry are actually collective exci- 
tations of atoms. This generically implies that their dispersion relation will 
be a relativistic one only at low energies (large scales), 21 and in each case 
there will be some short length scale (e.g., intermolecular distance for a fluid, 
coherence length for a superfluid, healing length for a BEC) beyond which 
deviations will be non-negligible. In general such microphysics induced cor- 
rections to the dispersion relation take the form 

E 2 = c 2 (m 2 c 2 + k 2 + A{k, K)) (247) 

where K is the scale that describes the transition to the full microscopic 
system (what we might call the "analogue Planck scale" ) . 

In general the best one can do is to expand A(k,K) around k = 0, so 
obtaining an infinite power series (of which it will be safe to retain only 
the lowest order terms), although in some special models (like BEC) the 
series is automatically finite due to intrinsic properties of the system. [In 
any case one can see that most of the analogue models so far considered 
lead to modifications of the form ±k 3 /K 2 or ±k A /K 2 .} Depending on the 
sign in front of the modification the group velocity at high energy can be 
larger (+) or smaller (— ) than the low energy speed of light c. These cases 
are usually referred in the literature as "superluminal" and "subluminal" 
dispersion relations. 

21 Actually, even relativistic behaviour at low energy can be non-generic, but we assume 
in this discussion that an analogue model by definition is a system for which all the 
linearized perturbations do propagate on the same Lorentzian background at low energies. 
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Most of the work on the trans-Planckian problem in the nineties focussed 
on studying the effect on Hawking radiation due to such modifications of 
the dispersion relations at high energies in the case of acoustic analogues 
188], |189| , |379| , |380| , |90| , and the question of whether such phenomenology 



could be applied to the case of real black holes (see e.g., [|^, |190| , |90| , |301|1 ). 22 
In all the aforementioned works Hawking radiation can be recovered under 
some suitable assumptions, as long as neither the black hole temperature nor 
the frequency at which the spectrum is considered are too close to the scale 
of microphysics K. However, the applicability of these assumptions to the 
real case of black hole evaporation is an open question. Also, in the case 
of the analogue models the mechanism by which the Hawking radiation is 
recovered is not always the same. We concisely summarize here the main 
results (but see e.g., ||382|| for further details). 



Subluminal dispersion relations: This was the case originally consid- 
ered by Unruh 1373), |5gQfl, 

u = K (tanh(k/K) n ) 1/n , (248) 

and by Corley and Jacobson |Xj 

to 2 = k 2 - k*/K 2 , (249) 

where both dispersion relations are given in the co-moving frame. 

The key feature is that in the presence of a subluminal modification the 
group velocity of the modes increases with k only up to some turning point 
(which is equivalent to saying that the group velocity does not asymptote to 
c, which could be the speed of sound, but instead is upper bounded). For 
values of k beyond the turning point the group velocity decreases to zero 
(for (|248|) ) or becomes imaginary (for (|249|) ) . In the latter case this can be 
interpreted as signifying the breakdown of the regime where the dispersion 
relation Q249|) can be trusted. The picture that emerged from these analy- 
ses concerning the origin of the outgoing Hawking modes at infinity is quite 
surprising. In fact, if one traces back in time an outgoing mode, as it ap- 
proaches the horizon it decreases its group velocity below the speed of sound. 



22 However see also [320, 324 for a radically different alternative approach based on 
the idea of "superoscillations" where ultrahigh frequency modes near the horizon can 
be mimicked (to arbitrary accuracy) by the exponential tail of an exponentially large 
amplitude mostly hidden behind the horizon. 
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At some point before reaching the horizon, the outgoing mode will appear 
as an ingoing one dragged into the black hole by the flow. Stepping further 
back in time it is seen that such a mode was located at larger and larger dis- 
tances from the horizon, and tends to a very high energy mode far away at 
early times. We have then that in this case one finds what might be called a 
"mode conversion" , where the origin of the outgoing Hawking quanta seems 
to originate from ingoing modes which have "bounced off" the horizon with- 
out reaching trans-Planckian frequencies in its vicinities. Several detailed 
analytical and numerical calculations have shown that such a conversion in- 
deed happens §80|, 0, 0, ||, |174|, |332|, [382| and that the Hawking result 
can be recovered for k <C K where k is the black hole surface gravity. 



Superluminal dispersion relations: The case of a superluminal disper- 
sion relation is quite different and, as we have seen, it has some experimental 
interest given that these are the kind of dispersion relations that arise in 
some promising analogue models (e.g., BECs). In this situation, the outgo- 
ing modes are actually seen as originating from behind the horizon. This 
implies that these modes somehow originate from the singularity (which can 
be a region of high turbulence in acoustic black hole analogues) , and hence 
it would seem that not much can be said in this case. However it is possible 
to show that if one imposes vacuum boundary conditions on these modes 
near the singularity, then it is still possible to recover the Hawking result, 
i.e., thermal radiation outside the hole fl89|. It is particularly interesting to 
note that this recovering of the standard result is not always guaranteed in 
the presence of superluminal dispersion relations. Corley and Jacobson |92| 
in fact discovered a very peculiar type of instability due to such superlumi- 
nal dispersion in the presence of black holes with inner horizons. The net 
result of the investigation carried out in |)2 is that the compact ergo-region 
characterizing such configurations is unstable to self-amplifying Hawking ra- 
diation. The presence of such an instability seems to be confirmed by the 
analysis carried on in [|138| , |139| where a Bose-Einstein condensate analogue 
black hole was considered. 



General conditions for Hawking radiation Is it possible to reduce the 
just described rather complex phenomenology to a few basic assumptions 
that must be satisfied in order to recover Hawking radiation in the pres- 
ence of Lorentz violating dispersion relations? A tentative answer is given 
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in [] 2J, where the robustness of the Hawking result is considered for gen- 
eral modified (subluminal as well as superluminal) dispersion relations. The 
authors of |[ • - /|| assume that the geometrical optics approximation breaks 
down only in the proximity of the event horizon (which is equivalent to say- 
ing that the particle production happens only in such a region). Here, the 
would be trans-Planckian modes are converted in sub-Planckian ones. Then, 
they try to identify the minimal set of assumptions that guarantees that 
such "converted modes" are generated in their ground states (with respect 
to a freely falling observer), as this is a well known condition in order to 
recover Hawking's result. They end up identifying three basic assumptions 
that guarantee such emergence of modes in the ground state at the horizon. 
First, the preferred frame selected by the breakdown of Lorentz invariance 
must be the freely falling one instead of the rest frame of the static observer 
at infinity (which coincides in this limit with the laboratory observer). Sec- 
ond, the Planckian excitations must start off in the ground state with respect 
to freely falling observers. Finally, they must evolve in an adiabatic way (i.e., 
the Planck dynamics must be much faster than the external sub-Planckian 
dynamics). Of course, although several systems can be found in which such 
conditions hold, it is also possible to show [|382j| that realistic situations in 
which at least one of these assumptions is violated can be imagined. It is 
hence still an open question whether real black hole physics does indeed sat- 
isfy such conditions, and it is hence robust against modifications induced by 
the violation of Lorentz invariance. 

Open issues In spite of the remarkable insight given by the above dis- 
cussed models (based on modified dispersion relations) it is not possible to 
consider them fully satisfactory in addressing the trans-Planckian problem. 



In particular it was soon recognized 91, 191] that in this framework it is not 



possible to explain the origin of the short wavelength incoming modes which 
are "progenitors" of the outgoing modes after bouncing off in the proximity 
of the horizon. For example, in the Unruh model (|248|) , one can see that if 



one keeps tracking an incoming mode back in time, then its group velocity 
(in the co-moving frame) drops to zero as its frequency gets more and more 
blue shifted. This is tantamount to saying that the trans-Planckian problem 
has just been moved from the region near the horizon out to the region near 
infinity. In the Corley-Jacobson model ( |. ; . | ) this unphysical behaviour is 
removed thanks to the presence of the physical cutoff K. However it is still 
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true that in tracking the incoming modes back in time one finally sees a wave 
packet so blue shifted that \k\ = K. At this point one can no longer trust the 
dispersion relation ( [ . | ) (which anyway in realistic analogue models is emer- 
gent and not fundamental), and hence the model has no predictive power 
regarding the ultimate origin of the relevant incoming modes. 

These conclusions regarding the impossibility of clearly predicting the 
origin at early times of the modes ultimately to be converted into Hawk- 
ing radiation is not specific to the particular dispersion relations ( | . | ) or 
(|2 19p one is using. The Killing frequency is in fact conserved on a static 
background, thus the incoming modes must have the same frequency as the 
outgoing ones, hence there can be no mode-mixing and particle creation. 
This is why one has actually to assume that the WKB approximation fails 
in the proximity of the horizon and that the modes are there in the vacuum 
state for the co-moving observer. In this sense the need for these assump- 
tions can be interpreted as evidence that these models are not fully capable 
of solving the trans-Planckian problem. 



Solid state and lattice models It was to overcome this type of issue that 
alternative ways of introducing an ultra-violet cutoff due to the microphysics 
were considered [ 20, c 21, 91]. In particular in [ the transparency of 



the refractive medium at high frequencies has been used to introduce an 
effective cutoff for the modes involved in Hawking radiation in a classical 
refractive index analogue model (see section |4.1.3| of this review). In this 
model an event horizon for the electromagnetic field modes can be simulated 
by a surface of singular electric and magnetic permeabilities. This would 
be enough to recover Hawking radiation but it would imply the unphysical 
assumption of a refractive index which is valid at any frequency. However 
it was shown in ||321|1 that the Hawking result can be recovered even in the 
case of a dispersive medium which becomes transparent above some fixed 
frequency K (which we can imagine as the plasma frequency of the medium), 
the only (crucial) assumption being again that the "trans-Planckian" modes 
with k > K are in their ground state near the horizon. 

There a lattice de- 
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An alternative avenue was instead considered in 
scription of the background was used for imposing a cutoff in a more physical 
way with respect to the continuum dispersive models previously considered. 
In such a discretized spacetime the field takes values only at the lattice points, 
and wavevectors are identified modulo 2n / 1 where I is the lattice character- 
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istic spacing, correspondingly one obtains a sinusoidal dispersion relation for 
the propagating modes. Hence the problem of recovering a smooth evolution 
of incoming modes to outgoing one is resolved by the intrinsically regularized 



behaviour of the wave vectors field. In |91] the authors explicitly considered 
the Hawking process for a discretized version of a scalar field, where the lat- 
tice is associated to the free-fall coordinate system (taken as the preferred 
system). With such a choice it possible to preserve a discrete lattice spacing. 
Furthermore the requirement of a fixed short distance cutoff leads to the 
choice of a lattice spacing constant at infinity, and that the lattice points are 
at rest at infinity and fall freely into the black hole. 23 In this case the lattice 
spacing grows in time and the lattice points spread in space as they fall to- 
ward the horizon. However this time dependence of the lattice points is found 
to be of order 1/k, and hence unnoticeable to long wavelength modes and rel- 
evant only for those with wavelengths of the order of the lattice spacing. The 
net result is that on such a lattice long wavelength outgoing modes are seen 
to originate from short wavelength incoming modes via a process analogous 



to the Bloch oscillations of accelerated electrons in crystals [OT, 191 



5.1.3 Horizon stability 

Although closely related, as we will soon see, we have to distinguish carefully 
between the mode analysis of a linear field theory (with or without modi- 
fied dispersion relations — MDR) over a fixed background and the stability 
analysis of the background itself. 

Let us consider a three-dimensional irrotational and inviscid fluid system 
with a stationary sink- type of flow (see figures p] and ||). The details of the 
configuration are not important for the following discussion, only the fact that 
there is a spherically symmetric fluid flow accelerating towards a central sink, 
that sink being surrounded by a sphere acting as a sonic horizon). Then, as 
we have discussed in section |^, linearizing the Euler and continuity equations 
leads to a massless scalar field theory over a black-hole-like spacetime. (We 
are assuming that the hydrodynamic regime remains valid up to arbitrarily 
short length scales; for instance, we are neglecting the existence of MDR). To 
be specific, let us choose the geometry of the canonical acoustic black-hole 



spacetime described in [ 



23 Referencc j9l| also considered the case of a lattice with proper distance spacing con- 
stant in time but this has the problem that the proper spacing of the lattice goes to zero 
at spatial infinity, and hence there is no fixed short distance cutoff. 
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ds 2 = -c 2 (l - dr 2 1 dr 2 + r 2 (d# 2 + sin 2 9 dp 2 ) . (250) 

In this expression we have used the Schwarzschild time coordinate r instead 
of the lab time t; c is constant. If we expand the field in spherical harmonics, 



<f>im(T,r,0,(p) = e 



Xim{r) 



r 



Y lm (0,<p), 



we obtain the following equation for the radial part of the field: 
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Here 

r* = r — (ro/4){ln(r + ro)/(r 

is a "tortoise" coordinate. 

In a normal mode analysis one requires boundary conditions such that 
the field is regular everywhere, even at infinity. However, if one is analyz- 
ing the solutions of the linear field theory as a way of probing the stability 
of the background configuration, one can consider less restrictive boundary 
conditions. For instance, one can consider the typical boundary conditions 
that lead to quasinormal modes: These modes have to be purely out-going 
at infinity and purely in-going at the horizon, but one does not require, for 
example, the modes to be normalizable. The quasinormal modes associated 
to this sink configuration have been analyzed in j31]|. The results found are 
qualitatively similar to those in the classical linear stability analysis of the 
Schwarzschild black hole in general relativity Q 55, 386 , |319| , | 52|, |269| . Of 
course, the gravitational field in general relativity has two dynamical de- 
grees of freedom — those associated with gravitational waves — that have to 
be analyzed separately; these are the "axial" and "polar" perturbations. In 
contrast, in the present situation we only have scalar perturbations. Nev- 
ertheless, the potentials associated with "axial" and "polar" perturbations 
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of Schwarzschild spacetime, and that associated with scalar perturbations of 
the canonical acoustic black hole, produce qualitatively the same behaviour: 
There is a series of damped quasinormal modes — proving the linear stability 
of the system — with higher and higher damping rates. 

An important point we have to highlight in here is that although in the 
linear regime the dynamical behaviour of the acoustic system is similar to 
general relativity, this is no longer true once one enters the non-linear regime. 
The underlying nonlinear equations in the two cases are very different. The 
differences are so profound, that in the general case of acoustic geometries 
constructed from compressible fluids, there exist sets of perturbations that, 
independently of how small they are initially, can lead to the development of 
shocks, a situation completely absent in vacuum general relativity. 

Now, given an approximately stationary, and at the very least metastable, 
classical black-hole-like configuration, a standard quantum mode analysis 
leads to the existence of Hawking radiation in the form of phonon emis- 
sion. This shows, among other things, that quantum corrections to the 
classical behaviour of the system must make the configuration with a sonic 
horizon dynamically unstable against Hawking emission. Moreover, in an 
analogue system with quantum fluctuations that maintain strict adherence 
to the equivalence principle (no MDR) it must then be impossible to create 
an isolated truly stationary horizon by external means — any truly station- 
ary horizon must be provided with an external power source to stabilize it 
against Hawking emission. That is, in an analogue system one could in princi- 
ple, by manipulating external forces, compensate for the backreaction effects 
that in a physical general relativity scenario cause the horizon to shrink (or 
evaporate) and thus becomes non-stationary. 

Let us describe what happens when one takes into account the existence 
of MDR. A wonderful physical system that has MDR explicitly incorporated 
in its description is the Bose-Einstein condensate. The macroscopic wave 
function of the BEC behaves as a classical irrotational fluid but with some 
deviations when short length scales become involved. (For length scales of 
the order of or shorter than the healing length.) What are the effects of the 
MDR on the dynamical stability of a black-hole like configuration in a BEC? 
The stability of a sink configuration in a BEC has been analyzed in [|138|, [139 



but taking the flow to be effectively one-dimensional. What they found is 
that these configurations are dynamically unstable: There are modes satis- 
fying the appropriate boundary conditions such that the imaginary parts of 
their associated frequencies are positive. These instabilities are associated 
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basically with the bound states inside the black hole. The dynamical ten- 
dency of the system to evolve is suggestively similar to that in the standard 
evaporation process of a black hole in semiclassical general relativity. This 
observation alone could make us question whether the first signatures of a 
quantum theory underlying general relativity might show up in precisely this 
manner. Interest in this question is reinforced by a specific analysis in the 
"loop quantum gravity" approach to quantizing gravity that points towards 



the the existence of fundamental MDR at high energies ||136|| . The formu- 
lation of effective gravitational theories that incorporate some sort of MDR 
at high energies is currently under investigation (see for example ||249| , ^]); 



these exciting developments are however beyond the scope of this review. 

Before continuing with the discussion on the stability of configurations 
with horizons, and in order not to cause confusion between the different 
wording used when talking about the physics of BECs and the emergent 
gravitational notions on them, let us write down a quite loose but useful 
translation dictionary: 

• The "classical" or macroscopic wave function of the BEC represents 
the classical spacetime of GR, but only when probed at long enough 
wavelengths. 

• The "classical" long-wavelength perturbations to a background solution 
of the Gross-Pitaevskii equation correspond to classical gravitational 
waves in GR. Of course this analogy does not imply that these are spin 
2 waves, it only points out that the perturbations are made from the 
same "substance" as the background configuration itself. 

• The macroscopic wave function of the BEC, without the restriction of 
being probed only at long wavelengths, corresponds to some sort of 
semiclassical vacuum gravity. Its "classical" behaviour (in the sense 
that does not involve any probability notion) is already taking into 
account its underlying quantum origin. 

• The Bogoliubov quantum quasiparticles over the "classical" wave func- 
tion correspond to a further step away from semiclassical gravity in 
that they are analogous to the existence of quantum gravitons over a 
(semiclassical) background spacetime. 

At this point we would like to remark, once again, that the analysis 
based on the evolution of a BEC has to be used with care. For example, they 
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cannot directly serve to shed light on what happens in the final stages of the 
evaporation of a black hole, as the BEC does not fulfil, at any regime, the 
Einstein equations. 

Now continuing the discussion, what happens when treating the pertur- 
bations to the background BEC flow as quantum excitations (Bogoliubov 
quasiparticles)?. What we certainly know is that the analysis of modes in a 
collapsing-to-form-a-black-hole background spacetime leads to the existence 



of radiation emission very much like Hawking emission [50, |93|], (this is 
why it is said that Hawking's process is robust against modifications of the 
physics at high energies). The comparison of these calculations with that 
of Hawking, (without MDR), tells us that the main modification to Hawk- 
ing's result is that now the Hawking flux of particles would not last forever 
but would vanish after a long enough time (this is why, in principle, we 
can dynamically create a configuration with a sonic horizon in a BEC). The 
emission of quantum particles reinforces the idea that the supersonic sink 
configurations are unstable. 
Summarizing: 

• If the perturbation to the BEC background configuration have "classi- 
cal seeds" (that is, are describable by the linearized Gross-Pitaevskii 
equation alone), then, one will have "classical" instabilities; 

• If the perturbations have "quantum seeds" (that is, are described by 
the Bogoliubov equations), then, one will have "quantum" instabilities. 

In the light of the acoustic analogies it is natural to ask whether there 
are other geometric configurations with horizons of interest, besides the sink 
type of configurations (these are the most similar to the standard description 
of black holes in general relativity, but probably not the simplest in terms 



of realizability in a real laboratory; for an entire catalogue of them see 0]). 
Here, we are going to specifically mention two effectively one-dimensional 
black hole-white hole configurations, one in a straight line and one in a ring 



(see figures [11], [12], respectively). 

A quantum mode analysis of the black hole-white hole configuration in 
a straight line, taking into account the existence of superluminal dispersion 
relations (similar to those in a BEC), led to the conclusion that the emission 
of particles in this configuration proceeds in a self-amplified (or runaway) 



manner [Q . We can understand this effect as follows: At the black horizon 
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Figure 11: One-dimensional velocity profile with a black hole horizon and a 
white hole horizon. 



a virtual pair of phonons are converted into real phonons, the positive en- 
ergy phonon goes towards infinity while the negative energy pair falls beyond 
the black horizon. However, the white horizon makes this negative energy 
phonon bounce back towards the black horizon (thanks to superluminal mo- 
tion) stimulating the emission of additional phonons. Although related to 
Hawking's process this phenomenon has a quite different nature: For ex- 
ample, there is no temperature associated with it. A stability analysis of a 
configuration like this in a BEC would lead to strong instabilities. This same 
configuration, but compactified into a ring configuration, has been dynam- 
ically analyzed in ||138| , |139|| . What they found is that there are regions of 



stability and instability depending on the parameters characterizing the con- 
figuration. We suspect that the stability regions appear because of specific 
periodic arrangements of the modes around the ring. Among other reasons, 
these arrangements are interesting because they could be easier to create 
in laboratory with current technology and their instabilities easier to detect 
than Hawking radiation itself. 

To conclude this subsection, we would like to highlight that there is still 
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Figure 12: One- dimensional velocity profile in a ring; the fluid flow exhibits 
two sonic horizons, one of black hole type and the other of white hole type. 

much to be learned by studying the different levels of description of an ana- 
logue system, and how they influence the stability or instability of configu- 
rations with horizons. 

5.2 Super-radiance 

Another phenomena that has been (and is being) analyzed from the analogue 
gravity perspective is super-radiance. The rotational kinetic energy accumu- 
lated in a rotating black hole can be extracted from it by scattering into it 
waves of sufficiently low frequency and high angular momentum. In general, 
in order that the wave can extract rotational energy from the system it has 
to satisfy the condition 



lo < m n 



(255) 
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where f2 is the angular speed of the black hole at the event horizon and 
m is the harmonic azimuthal number of the wave. This is a purely classi- 



cal process first considered by Penrose [307]. When dealing with quantum 



fields, as opposed to classical fields, this process can proceed spontaneously: 
Quantum mechanically, a rotating black hole will tend to radiate away all of 
its angular momentum, eventually approaching a non-rotating Schwarzschild 
black hole |[430| , [431|| . This process is known as super-radiance. (The term 
super-radiance was already used in condensed matter to describe processes 
in which there was some coherent emission of radiation from an otherwise 
disordered system). 

Again, these processes have a purely kinematical origin so they are per- 
fectly suitable for being reproduced in an analogue model. Regarding these 
processes, the simplest geometry that one can reproduce, thinking of ana- 
logue models based on fluid flows, is that of the draining bathtub of Section 
0. Of course, this metric does not exactly correspond to Kerr geometry, nor 



even to a section of it ||402|| . However, it is qualitatively similar. It can be used 
to simulate both Penrose's classical process and quantum super-radiance as 
these effects do not depend on the specific multipole decomposition of Kerr's 
geometry, but only on its rotating character. A specific experimental setup 
has been put forward by Schutzhold and Unruh using gravity waves in a shal- 
low basin mimicking an ideal draining bathtub |347|| . Equivalently to what 
happens with Kerr black holes, this configuration is classically stable (in the 
linear regime) | T|. A word of caution is in order here: Interactions of the 
gravity surface waves with bulk waves (neglected in the analysis) could cause 
the system to became unstable ||417]| . This instability has no counterpart in 
standard general relativity (though it might have one in braneworld theo- 
ries) . Super-resonant scattering of waves in this rotating sink configuration, 
or in a simple purely rotating vortex, could in principle be observed in this 
and other analogue models. There are already several articles dealing with 



this problem pi, [27], |26|, |64|, [TT4] p9 



A related phenomenon one can consider is the black-hole bomb mech- 
anism ||314|| . One would only have to surround the rotating configuration 
by a mirror for it to become grossly unstable. What causes the instability 
is that those in-going waves that are amplified when reflected in the ergo- 
sphere would then in turn be reflected back toward the ergoregion, due to 
the exterior mirror, thus being amplified again, and so on. 

An interesting phenomenon that appears in many condensed matter sys- 
tems is the existence of quantized vortices. The angular momentum of these 
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vortices come in multiples of some fundamental unit (typically H or some- 
thing proportional to H). The extraction of rotational energy by Penrose 
process in these cases could only proceed via finite-energy transitions. This 
would supply an additional specific signature to the process. In such a highly 
quantum configuration it is also important to look for the effect of having 
high-energy dispersion relations. For example, in BECs, the radius of the 
ergoregion of a single quantized vortex is of the order of the healing length, 
so one cannot directly associate an effective Lorentzian geometry to this por- 
tion of the configuration. Any analysis that neglects the high energy terms 
is not going to give any sensible result in these cases. 



5.3 Cosmological geometries 

Analogue model techniques have also been applied to cosmology. In a cosmo- 
logical framework the key items of interest are the Friedmann-Robertson- 
Walker [FRW] geometries, more properly called the Friedmann-Lemaitre- 
Robert son- Walker [FLRW] geometries. The simulation of such geometries 
has been considered in various works such as [|17], pl| |117| , |116| , |115|, [244| , [58| , 



9[ |42 5|| with a specific view to enhancing our understanding of "cosmological 
particle production" driven by the expansion of the universe. 
The acoustic metric can be written as 



ds 1 



dt 2 -2v-dxdt + dx 2 ] 



(256) 



Essentially there are two ways to use this metric to reproduce cosmological 
spacetimes: One is based on physical explosion, the other on rapid variations 
in the "effective speed of light" . 

Following fI8| , |116| , |205| , |206|| one can take an homogeneous system 
p(t), c s (t) and a radial profile for the velocity v = (b/b)r, with b a scale factor 
depending only on t. Then, defining a new radial coordinate as rj = r/6 the 
metric can be expressed as 



As' 



dt 2 



+ r 2 b dQj 



(257) 



Introducing a Hubble-like parameter, 



H b (t) 



b( t y 



(258) 
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the equation of continuity can be written as 



p + 3H b (t) p = 0; 



Pit) 



Po 



(259) 



with po constant. Finally we arrive at the metric of a flat FLRW geometry 



with 



ds 2 = _ T 2 (t) dt 2 + fl 2 (t) (dr 2 + r 2 dQ 2^ ^ 



(261) 



T(t) = y/pc~ a \ a s (t) = J—b. 

V 

The proper Friedmann time, r, is related to the laboratory time, t, by 



T 



T(t) dt. 



The other avenue starts from a fluid at rest v 
laboratory at all times: 



(262) 



with respect to the 



ds 2 



-pc s dt z 



^dx 2 . 



(263) 



Now it is not difficult to imagine situation in which p remains constant, in 
a sufficiently large region of space, while the speed of sound decreases with 
time (this can be made in BECs for example by changing with time the 
value of the scattering length |L7|, |I8|]). This again reproduces a expanding 
flat FLRW Universe. 

Models considered to date focus on variants of the BEC-inspired ana- 
logues: 



Fedichev and Fischer [ 116 , 1 [] have investigated WKB estimates of 
the cosmological particle production rate and (1 + 1) dimensional cos- 
mologies, both in expanding BECs. 



Lidsey p44|| , and Fedichev and Fischer [|117|| have focussed on the be- 
haviour of cigar-like condensates in grossly asymmetric traps. 

Barcelo et al. [17], [l^] have focussed on the central region of the BEC 
and thereby tried (at least locally) to mimic FLRW behaviour as closely 
as possible. 
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In [ | Fischer and Schutzhold propose the use of two-component 
BECs to simulate cosmic inflation. 



Weinfurtner [ |26 > [425| has concentrated on the approximate simulation 
of de Sitter spacetimes. 



In all of these models the general expectations of the relativity community 
have been borne out — theory definitely predicts particle production, and the 
very interesting question is the extent to which the formal predictions are 
going to be modified when working with real systems experimentally fl8| . 
We expect that these analogue models provide us with new insights as to 
how their inherent modified dispersion relations affect cosmological process 
such as the generation of a primordial spectrum of perturbations (see for 
example |§ |3§ g§, gt| g7[ [71], |TUB|, [TBD], git], f^T} |232|, [25T], ^ [ZS3|, 

|2544 \276[ \277[ |351| , |363|| where analogue-like ideas are applied to cosmological 
inflation). 

An interesting side-effect of the original investigation, is that birefringence 
can now be used to model "variable speed of light" [VSL] geometries |109 



Since analogue models quite often lead to two or more "excitation cones", 
rather than one, it is quite easy to obtain a bimetric or multi-metric model. If 
one of these metric is interpreted as the "gravitational" metric and the other 
as the "photon" metric, then VSL cosmologies can be given a mathematically 
well-defined and precise meaning PR |109 . 



5.4 Bose novae: an example of the reverse flow of in- 
formation? 

As we have seen in the previous sections, analogue models have in the past 
been very useful in providing new, condensed matter physics inspired, ideas 
about how to solve longstanding problems of semiclassical gravity. In closing 
this section it is interesting to briefly discuss what perhaps represents, so far, 
the only tentative attempt to use analogue models in the reverse direction; 
that is to import well known concepts of semiclassical gravity into condensed 
matter frameworks. 

The phenomenon we are referring to is the so called "Bose nova" |105 . 
This is an experiment dealing with a gas of a few million 85 Rb atoms at a 
temperature of about 3 nK. The condensate is rendered unstable by exploit- 
ing the possibility of tuning the interaction (more precisely the scattering 
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length) between the atoms via a magnetic field. Reversing the sign of the 
interaction, making it attractive, destabilizes the condensate. After a brief 
waiting time (generally called tcoiiapse), the condensate implodes and loses a 
sizeable fraction of its atoms in the form of a "nova burst" . If left to evolve 
undisturbed, the number of atoms in the burst stabilizes and a remnant con- 
densate is left. However if the condensate interaction is again made repulsive 
after some time t e voive, before the condensate has sufficient time to stabilize, 
then the formation of "jets" of atoms is observed, these jets being character- 
ized by lower kinetic energy and a distinct shape with respect to the burst 
emission. 

Interestingly, an elegant explanation of such a phenomenology was pro- 
posed in [ p8| , |59[] , based on the well known semiclassical gravity analysis of 
particle creation in an expanding universe. In fact the dynamics of quan- 
tum excitations over the collapsing BEC was shown to closely mimic that 
for quantum excitations in a time-reversed (collapsing instead of expanding) 
scenario for cosmological particle creation. This is not so surprising as the 
quantum excitations above the BEC ground state feel a time-varying back- 
ground during the collapse, and as a consequence one then expects squeezing 
of the vacuum state and mode mixing which are characteristic of quantum 
field theory in variable external fields. 

However the analogy is even deeper than this. In fact in ]5^, [59| a key 



role in explaining the observed burst and jets is played by the cosmology 
borrowed concepts of "frozen" versus "oscillating" modes (although with a 
reverse dynamics with respect to the standard [expanding] cosmological case). 
In the case of Bose novae the modes which are amplified are those for which 
the physical frequency is smaller than the collapse rate, while modes with 
higher frequencies remain basically unaffected and their amplitudes obey a 
harmonic oscillator equation. As the collapse rate decreases, more and more 
modes stop growing and start oscillating, which is equivalent to a creation 
of particles from the quantum vacuum. In the case of a sudden stop of the 
collapse by a new reversal of the sign of the interaction all of the previously 
growing modes are suddenly converted into particles, explaining in this way 
the generation of jets and their lower energy (they correspond to modes with 
lower frequencies with respect to those generating the bursts). 

Although this simple model cannot explain all the details of the Bose 
novae phenomenology, we think it is remarkable how far it can go in explain- 
ing several observed features by exploiting the language and techniques so 
familiar to quantum cosmology. In this sense the analysis presented in 5J| 
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primarily shows a possible new application of analogue models where they 
could be used to lend ideas and techniques developed in the context of grav- 
itational physics to the explanation of condensed matter phenomena. 



5.5 Going further 



For more details on the trans-Planckian problem some of the key papers 
are the relatively early papers of Unruh ||379| , |380|| , and Jacobson ||188| , |189 



For super-radiance and cosmological issues there seems to be considerable 
ongoing interest, and one should carefully check Spires for the most recent 
articles. 
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6 Future directions 



The key question one should ask at this stage is: "Where can we go from 
here?" Apart from continuing with the analysis of the issues described in 
the previous section, there are additional interesting routes, both theoretical 
and experimental, worth of future exploration. In particular, the following 
topics come to mind: 

• Back reaction, 

• Equivalence principle, 

• Emergent gravity, 

• Quantum gravity phenomenology, 

• Quantum gravity, 

• Experimental analogue gravity. 

There has been already some works dealing with this topics in the contest of 
analogue gravity. Let us now expand on these issues a little. 



6.1 Back reaction 

There are important phenomena on gravitational physics whose understand- 
ing needs of analysis well beyond classical general relativity and field theory 
on curved background spacetimes. The black hole evaporation process can 
be considered as paradigmatic among these phenomena. Here, we partic- 
ularize our discussion to this case. Since we are currently unable to ana- 
lyze the entire process of black hole evaporation within a complete quantum 
theory of gravity, a way of proceeding is to analyze the simpler (but still 
extremely difficult) problem of semiclassical back reaction (see for exam- 
ple ||101| , |80| , |34[ , |134| , |5T| , |259|| ). One takes a background black hole spacetime, 



calculates the expectation value of the quantum energy-momentum tensor of 
matter fields in the appropriate quantum state (the Unruh vacuum state for 
a radiating black hole), and then takes this expectation value as a source 
for the perturbed Einstein equations. This calculation gives us information 
about the tendency of spacetime to evolve under vacuum polarization effects. 

A nice feature of analogue models of general relativity is that although 
the underlying classical equations of motion have nothing to do with Einstein 
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equations, the tendency of the analogue geometry to evolve due to quantum 
effects is formally equivalent (approximately, of course) to that in semiclassi- 
cal general relativity. Therefore, the onset of the back reaction effects (if not 
their precise details) can be simulated within the class of analogue models. 
An example of the type of back-reaction calculations one can perform are 
those in ||, [TfJ. These authors started from an effectively one-dimensional 
acoustic analogue model, configured to have an acoustic horizon by using a 
Laval nozzle to control the flow's speed. They then considered the effect of 
quantizing the acoustic waves over the background flow. To calculate the ap- 
propriate back reaction terms they took advantage of the classical conformal 
invariance of the (l+l)-dimensional reduction of the system. In this case, we 
know explicitly the form of the expectation value of the energy-momentum 
tensor trace (via the trace anomaly). The other two independent components 
of the energy-momentum tensor were approximated by the Polyakov stress 
tensor. In this way, what they found is that the tendency of a left-moving 
flow with one horizon is for it to evolve in such a manner as to push the hori- 
zon down-stream at the same time that its surface gravity is decreased. This 
is a behaviour similar to what is found for near-extremal Reissner-Nordstrom 
black holes. (However, we should not conclude that acoustic black holes are 
in general closely related to near-extremal Reissner-Nordstrom black holes, 
rather than to Schwarzschild black holes. This result is quite specific to the 
particular one-dimensional configuration analyzed.) 

Can we expect to learn something new about gravitational physics by 
analysing the problem of back reaction in different analogue models? As we 
have repeatedly commented, the analyses based on analogue models force 
us to consider the effects of modified high-energy dispersion relations. For 
example, in BECs, they affect the "classical" behaviour of the background 
geometry as much as the behaviour of the quantum fields living on the back- 
ground. In seeking a semiclassical description for the evolution of the geome- 
try, one would have to compare the effects caused by the modified dispersion 
relations to those caused by pure semiclassical back reaction (which incor- 
porates deviations from standard general relativity as well). In other words, 
one would have to understand the differences between the standard back re- 
action scheme in general relativity, and that based on equations ( |18tj ) and 

13). 



To end this subsection we would like to comment that one can go beyond 
the semiclassical back-reaction scheme by using the so-called stochastic semi- 
classical gravity programme ||183| , |184| , |185|| . This programme aims to pave 



108 



the way from semiclassical gravity toward a complete quantum-gravitational 
description of gravitational phenomena. This stochastic gravity approach not 
only considers the expectation value of the energy-momentum tensor but also 
its fluctuations, encoded in the semiclassical Einstein-Langevin equation. In 
a very interesting paper ||301|| , Parentani showed that the effects of the fluc- 



tuations of the metric (due to the in-going flux of energy at the horizon) on 
the out-going radiation led to a description of Hawking radiation similar to 
that obtained with analogue models. It would be interesting to develop the 
equivalent formalism for quantum analogue models and to investigate the 
different emerging approximate regimes. 

6.2 Equivalence principle 

Analogue models are of particular interest in that the analogue spacetimes 
that emerge often violate, to some extent or another, the Einstein equivalence 
principle [H>, 400 1. Since the Einstein equivalence principle (or more precisely 



the universality of free fall) is experimentally tested to the accuracy of about 
1 part in 10 13 , it is important to build this principle into realistic models of 
gravitation — most likely at a fundamental level. 

One way of interpreting the Einstein equivalence principle is as a "princi- 
ple of universality" for the geometrical structure of spacetime. Whatever the 
spacetime geometrical structure is, if all excitations "see" the same geometry 
one is well on the way to satisfying the observational and experimental con- 
straints. In a metric theory, this amounts to the demand of mono-metricity: 
a single universal metric must govern the propagation of all excitations. 

Now it is this feature that is relatively difficult to arrange in analogue 
models. If one is dealing with a single degree of freedom then mono-metricity 
is no great constraint. But with multiple degrees of freedom, analogue space- 
times generally lead to refringence — that is the occurrence of Fresnel equa- 
tions that often imply multiple propagation speeds for distinct normal modes. 
To even obtain a bi-metric model (or more generally, a multi-metric model), 
requires an algebraic constraint on the Fresnel equation that it completely 
factorize into product of quadratics in frequency and wavenumber. Only 
if this algebraic constraint is satisfied can one assign a "metric" to each of 
the quadratic factors. If one further wishes to impose mono-metricity then 
the Fresnel equation must be some integer power of some single quadratic 



expression, an even stronger algebraic statement [16, 400 . 
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Faced with this situation, there are two ways in which the analogue grav- 
ity community might proceed: 

• Try to find a broad class of analogue models (either physically based or 
mathematically idealized) that naturally lead to mono-metricity. Little 
work along these lines has yet been done; at least partially because it is 
not clear what features such a model should have in order to be "clean" 
and "compelling". 

• Accept refringence as a common feature of the analogue models and 
attempt to use refringence to ones benefit in one or more ways: 

— There are real physical phenomena in non-gravitational settings 
that definitely do exhibit refringence and sometimes multi-metricity. 
Though situations of this type are not directly relevant to the grav- 
ity community, there is significant hope that the mathematical and 
geometrical tools used by the general relativity community might 
in these situations throw light on other branches of physics. 

— Use the refringence that occurs in many analogue models as a 
way of "breaking" the Einstein equivalence principle, and indeed 
as a way of "breaking" even more fundamental symmetries and 
features of standard general relativity, with a view to exploring 
possible extensions of general relativity. While the analogue mod- 
els are not themselves primary physics, they can nevertheless be 
used as a way of providing hints as to how more fundamental 
physics might work. 



6.3 Emergent gravity 

One of the more fascinating approaches to "quantum gravity" is the sug- 
gestion, typically attributed to Sakharov ||334| , |395|| , that gravity itself may 
not be "fundamental physics". Indeed it is now a relatively common opin- 
ion, maybe not mainstream but definitely a strong minority opinion, that 
gravity (and in particular the whole notion of spacetime and spacetime ge- 
ometry) might be no more "fundamental" than is fluid dynamics. The word 
"fundamental" is here used in a rather technical sense — fluid mechanics 
is not fundamental because there is a known underlying microphysics; that 
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of molecular dynamics, of which fluid mechanics is only the low-energy low- 
momentum limit. Indeed the very concepts of density and velocity field, 
which are so central to the Euler and continuity equations, make no sense at 
the microphysical level and emerge only as one averages over timescales and 
distance-scales larger than the mean free time and mean free path. 

In the same way, it is plausible (even though no specific and compelling 
model of the relevant microphysics has yet emerged) that the spacetime man- 
ifold and spacetime metric might arise only once one averages over suitable 
microphysical degrees of freedom. Sakharov had in mind a specific model in 
which gravity could be viewed as an "elasticity" of the spacetime medium, 



and was "induced" via one- loop physics in the matter sector [j 54], |395|| . In 
this way Sakharov had hoped to relate the observed value of Newton's con- 
stant (and the cosmological constant) to the spectrum of particle masses. 

More generally the phrase "emergent gravity" is now used to describe the 
whole class of theories in which the spacetime metric arises as a low-energy 
approximation, and in which the microphysical degrees of freedom might 
be radically different. Analogue models, and in particular analogue models 
based on fluid mechanics or the fluid dynamic approximation to BECs, are 
specific examples of "emergent physics" in which the microphysics is well 
understood. As such they are useful for providing hints as to how such a 
procedure might work in a more fundamental theory of quantum gravity. 

6.4 Quantum gravity — phenomenology 

Over the last few years a widespread consensus has emerged that observa- 
tional tests of quantum gravity are for the foreseeable future likely to be 
limited to precision tests of dispersion relations and their possible deviations 
from Lorentz invariance |[263|| . The key point is that at low energies (well 



below the Planck energy) one expects the locally Minkowskian structure of 
the spacetime manifold to guarantee that one sees only special relativistic 
effects; general relativistic effects are negligible at short distances. How- 
ever as ultra high energies are approached (although still below Planck scale 
energies) several quantum gravity models seem to predict that the locally 
Euclidean geometry of the spacetime manifold will break down. There are 
several scenarios for the origin of this breakdown ranging from string the- 



ory I \J\, |110|| to brane worlds | | and loop quantum gravity ||136|| . Common 
to all such scenarios is that the microscopic structure of spacetime is likely to 
show up in the form of a violation of Lorentz invariance leading to modified 
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dispersion relations for elementary particles. Such dispersion relations are 
characterized by extra terms (with respect to the standard relativistic form) 
which are generally expected to be suppressed by powers of the Planck en- 
ergy. Remarkably the last years have seen a large wealth of work in testing 
the effects of such dispersion relations and in particular strong constraints 
have been cast by making use of high energy astrophysics observations (see 
e.g. i, H, [197], [l96|, [198], [199], g6j| and references therein). 



Several of the analogue models are known to exhibit similar behaviour, 
with a low-momentum effective Lorentz invariance eventually breaking down 
at high momentum once the microphysics is explored. 24 Thus some of the 
analogue models provide controlled theoretical laboratories in which at least 
some forms of subtle high-momentum breakdown of Lorentz invariance can 
be explored. As such the analogue models provide us with hints as to what 
sort of modified dispersion relation might be natural to expect given some 
general characteristics of the microscopic physics. Hopefully investigation of 
appropriate analogue models might be able to illuminate possible mechanisms 
leading to this kind of quantum gravity phenomenology, and so might be able 
to provide us new ideas about other effects of physical quantum gravity that 
might be observable at sub-Planckian energies. 

6.5 Quantum gravity — fundamental models 

When it comes to dealing with "fundamental" theories of quantum gravity, 
the analogue models play an interesting role that is complementary to the 
more standard approaches. The search for a quantum theory of gravity is 
fundamentally a search for an appropriate mathematical structure in which to 
simultaneously phrase both quantum questions and gravitational questions. 
More precisely, one is searching for a mathematical framework in which to 
develop an abstract quantum theory that then itself encompasses classical 
Einstein gravity (the general relativity), and reduces to it in an appropriate 
limit p5j g58j plj. 



The three main approaches to quantum gravity currently in vogue, "string 
models" [also known as "M-models"], "loop space" [and the related "spin 
foams"], and "lattice models" [Euclidean or Lorentzian] all share one feature: 

It is however important to keep in mind that not all the above cited quantum gravity 
models violate the Lorentz symmetry in the same manner. The discreteness of spacctime 
at short scales is not the only way of breaking Lorentz invariance. 
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They attempt to develop a "pre-geometry" as a replacement for classical dif- 
ferential geometry (which is the natural and very successful mathematical 
language used to describe Einstein gravity) |H| |358| , |131| , |329| , |328| , ^0[ . The 



basic idea is that the mooted replacement for differential geometry would 
be relevant at extremely small distances (where the quantum aspects of the 
theory would be expected to dominate), while at larger distances (where the 
classical aspects dominate) one would hope to recover both ordinary differen- 
tial geometry and specifically Einstein gravity or possibly some generalization 
of it. The "string", "loop", and "lattice" approaches to quantum gravity dif- 
fer in detail in that they emphasise different features of the long-distance 
model, and so obtain rather different short- distance replacements for clas- 
sical differential geometry. Because the relevant mathematics is extremely 
difficult, and by and large not particularly well understood, it is far from 
clear which if any of these three standard approaches will be preferable in 
the long run ||358| . 



We feel it likely that analogue models can shed new light on this very con- 
fusing field by providing a concrete specific situation in which the transition 
from the short- distance "discrete" or "quantum" theory to the long-distance 
"continuum" theory is both well understood and non-controversial. Here 
we are specifically referring to fluid mechanics, where at short distances the 
system must be treated using discrete atoms or molecules as the basic build- 
ing blocks, while at large distances there is a well-defined continuum limit 
that leads to the Euler and continuity equations. Furthermore once one is 
in the continuum limit, there is a well-defined manner in which a notion of 
"Lorentzian differential geometry" , and in particular a "Lorentzian effective 



spacetime metric" can be assigned to any particular fluid flow [378, 391, 286 



Indeed, the "analogue gravity programme" is extremely successful in this re- 
gard, providing a specific and explicit example of a "discrete" — > "continuum" 
— > "differential geometry" chain of development. What the "analogue grav- 
ity programme" does not seem to do as easily is to provide a natural direct 
route to the Einstein equations of general relativity, but that merely indicates 
that current analogies have their limits and therefore, one should not take 
them too literally [378,283]. Fluid mechanics is a guide to the mathematical 
possibilities, not an end in itself. The parts of the analogy that do work well 
are precisely the steps where the standard approaches to quantum gravity 
have the most difficulty, and so it would seem useful to develop an abstract 
mathematical theory of the "discrete" — > "continuum" — > "differential ge- 
ometry" chain using this fluid mechanical analogy (and related analogies) as 
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inspiration. 

6.6 Going further 

Beyond the various theoretical issues we have discussed above there is the 
important question of "experimental analogue gravity" - - to what extent 
can all these ideas be tested in direct laboratory experiments? Currently 
several experimental groups are interested, but to the best of our knowledge 
no actual experiments are currently underway. Broadly speaking, for any 
experimental group interested in analogue spacetimes the two key issues to 
address are: 

• Identify a particular analogue model easily amenable to laboratory in- 
vestigation, and double check the extent to which the model provides 
a theoretically robust and clean analogue to general relativistic curved 
spacetime. 

• Identify the technical issues involved in actually setting up a laboratory 
experiment. 

While the consensus in the theoretical community is that Bose-Einstein 
condensates are likely to provide the best working model for analogue gravity, 
it is possible that we might still be surprised by experimental developments. 
We leave this as an open challenge to the experimental community. 
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7 Conclusions 



In this review article we have seen the interplay between standard general 
relativity and various analogies that can be used to capture aspects of its be- 
haviour. These analogies have ranged from rather general but very physical 
analogue models based on fluid-acoustics, geometrical optics, and wave op- 
tics, to highly specific models based on BECs, liquid helium, slow light, etc. 
Additionally, we have seen several rather abstract mathematical toy mod- 
els that bring us to such exotic structures and ideas such as birefringence, 
bimetricity, Finsler spaces, and Sakharov's induced gravity. 

The primary reason that these analogies were developed within the gen- 
eral relativity community was to help in the understanding of general relativ- 
ity by providing very down-to-earth models of otherwise subtle behaviour in 
general relativity. Secondary reasons include the rather speculative sugges- 
tion that there may be more going on than just analogy — it is conceivable 
[though perhaps unlikely] that one or more of these analogue models could 
suggest a relatively simple and useful way of quantizing gravity that side- 
steps much of the technical machinery currently employed in such efforts. 
A tertiary concern [at least as far as the general relativity community is 
concerned] is the use of relativity and differential geometric techniques to 
improve understanding of various aspects of condensed matter physics. 

The authors expect interest in analogue models to continue unabated, 
and suspect that there are several key but unexpected issues whose resolution 
would be greatly aided by the analysis of appropriate analogue models 

7.1 Going further 

Though every practicing scientist already know this, for the sake of any 
student reading this we mention the following resources: 

• http://www.slac.stanford.edu/spires/ — the bibliographic database for 
keeping track of (almost all of) the relevant literature. 

• http://www.arXiv.org — the electronic-preprint (e-print) database for 
accessing the text of (almost all, post 1992) relevant articles. 

• http: / /www. Iivingreviews.org/ — the Living Reviews portal. 

Those three access points should allow you to keep abreast of what is going 
on in the field. 
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